o 

(N 
X> 

PL, 

m 



-f— > 

CO 

a 

C 

o 

> 

(N 
(N 
(N 

(N 
O 



Cofermion Theory for Pseudogap Phenomena and Superconducting Mechanism of 

Underdoped Cuprate Superconductors 

Youhei YamajQ and Masatoshi Imada 
Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan 

(Dated: January 12, 2013) 

We study pseudogap phenomena and Fermi-arc formation experimentally observed in typical two 
dimensional doped Mott insulators, namely, underdoped cuprate superconductors. To develop a 
physically unequivocal theory, we start from the slave-boson mean-field theory for the Hubbard 
model on a square lattice. Our crucial step is to further take into account the charge dynamics 
and fluctuations. The extra charge fluctuations seriously modify low-energy single-particle spectra 
of doped Mott insulators near the Fermi level: An electron added around an empty site (or a hole 
added around a doubly occupied site) constitutes composite fermion (cofermion) , called holo-electron 
(or doublo-hole) at low energy in distinction from the normal quasiparticles. These unexplored 
composite fermions substantiate the extra charge fluctuation. We show that the quasiparticles 
hybridize with the holo-electrons and doublo-holes. The resultant hybridization gap is identified as 
the pseudogap observed in the underdoped region of the high-T c cuprates. Because the Fermi level 
crosses the top (bottom) of the low-energy band formed just below (above) the hybridization gap in 
the hole-doped (electron-doped) case, it causes a Fermi-surface reconstruction, namely, a topological 
change in the Fermi surface forced by the penetration of zeros of the quasiparticle Green function. 
This reconstruction signals the emergence of a non-Fermi-liquid phase. The pseudogap, and the 
resultant formation of pocket or arc of the Fermi surface reproduce the experimental results for the 
cuprate superconductors in the underdoped region. The pairing channel opens not only between 
two quasiparticles but also between a quasiparticle and a cofermion. This pairing solves the puzzle 
of the dichotomy between the d-wave superconductivity and the precursors of the the insulating gap 
in the antinodal region. We propose and analyze them as the mechanism of the high-temperature 
superconductivity for the cuprates. 



I. INTRODUCTION 

The discovery of cuprate superconductors has triggered 
extensive studies on the nature of low-energy electronic 
excitations evolving in doped Mott insulators. The ex- 
tensive interest on the doped Mott insulators exists be- 
cause it must be directly related with the origin of the 
high temperature superconductivity itself. 

From the early stage of the studies, experimental obser- 
vations have indicated that quasiparticle states of normal 
phases of cuprates change qualitatively with the increase 
of the doping, together with changes in superconducting 
transition temperatures. It has been showr>L2 that the 
Hall coefficient Rjj changes its sign near the so-called op- 
timal doping for the highest superconducting transition 
temperatures, followed by a steep increase of its ampli- 
tude \Rh\ with lowering doping concentration, typically 
indicating a drastic change in low-energy quasiparticles of 
the normal state between the underdoped and overdoped 
regions. Such a drastic change is also but differently 
suggested from observations of pseudogap phenomena^—, 
where the spin and charge excitations are unexpectedly 
suppressed in the underdoped region. Various types of 
non-Fermi-liquid properties are accompanied in this re- 
gion. 

Recent improved experimental tools have enabled re- 
solving low-energy single-particle spectra near the Fermi 
level. In particular, strongly momentum-dependent 
quasiparticle states in the hole-underdopcd cuprates2»i£ 
observed by angle-resolved photocmission spectroscopy 



(ARPES) studies have renewed the interest in the low- 
energy spectrum of the cuprate superconductors. In con- 
trast to the overdoped region, where a large Fermi surface 
crossing the region around the so-called antinodal points 
(±7r, 0) and (0, ±7r) in the 2D Brillouin zone for the Cu02 
plane is clearly observed, low-energy quasiparticle states 
around the antinodal points are missing in the under- 
doped cuprates. It emerges as a truncation of the large 
Fermi surface observed in the overdoped region. The re- 
sultant truncated structure is called the "Fermi arc" . 

LHB UHB 
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FIG. 1: Schematic global density of states of hole-doped 
Mott insulators. The coherent band around the Fermi level 
/U, the upper Hubbard band (UHB), and the lower Hubbard 
band (LHB) are schematically illustrated. The low-energy 
unoccupied spectral weight (LUSW) is also indicated. 



From more fundamental point of view, the normal 
state of the cuprates offers a challenge of condensed mat- 
ter physics as a typical open issue of "Mott physics", 
namely, nature of strongly correlated metals in the prox- 
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imity to the Mott insulator— ~—. The experimentally ob- 
served arc-like Fermi surface in the underdoped cupratcs 
is a hallmark of the proximity to the Mott insulators. 
As we illustrate in FigfTJ global energy spectra of doped 
Mott insulators are known to consist of essentially three 
energy "bands" ; a coherent band around the Fermi level 
(i and two incoherent bands, namely, the upper Hub- 
bard band (UHB) located above /j,, and the lower Hub- 
bard band (LHB) located below \i. For the hole-doped 
(electron-doped) systems, the coherent band is formed 
around the top of the LHB (the bottom of the UHB). 
The spectral weight formed just above (below) /i within 
the coherent band is often called the low-energy unoccu- 
pied spectral weight (LUSW) in the hole-doped (electron- 
doped) Mott insulators. 

Such a global structure of spectra can be roughly de- 
scribed by a simple picture given by the dynamical mean- 
field theory^. It has unified two scenarios of earlier stud- 
ies by HubbarcU^, and Brinkman and Ricei^ based on 
the Gutzwillcr approximation-^. The Hubbard approxi- 
mation captures the formation of UHB and LHB. In the 
Brinkman-Rice picture, Mott insulators appear as a con- 
sequence of homogeneous vanishing of the quasiparticle 
weight uniformly on the Fermi surface. Therefore, this 
approximation draws a picture that quasiparticles are 
renormalizcd in a momentum-independent fashion and 
the effective mass diverges independently of the momen- 
tum position on the Fermi surface on the verge of the 
Mott transition. 

In infinite-dimensional systems, such a route of the 
Mott transition has turned out to be correc t 14 ' 18 ' 19 . How- 
ever, in finite dimensions, this picture has turned out to 
be too simple to understand the low-energy spectra and 
the nature of the Mott transition: Realistic theory has to 
capture significant momentum-dependent quasiparticle 
renormalization, as has been experimentally observed as 
Fermi-arc formation in the underdoped cupratcs, and has 
also been revealed from various numerical attempts^— . 
It has now become increasingly clear that a conceptu- 
ally new idea for describing Mott physics including an 
emergence of the new non-Fcrmi-liquid phase in the un- 
derdoped region is desired. 

In this paper, we extend our previous attempt^ for 
describing the momentum-dependent Mott physics in- 
cluding the mechanism of Mott transitions at the level 
of LUSW and for understanding the experimentally sup- 
ported dramatic change in the electronic states of the 
cuprates. We give a detailed comparison with experi- 
ments and a novel superconducting mechanism based on 
the previous attempt 2 ^, together with the derivation of 
our theoretical description in details. We particularly 
pay attention to the possible reconstruction of the Fermi 
surface. If the reconstruction emerges, the doped Mott 
insulator becomes topologically inequivalent to the con- 
ventional Fermi liquid and hence offers an unexplored 
avenue for long-sort issue of unconventional metals. 

In the present work, we propose a simple and physi- 
cally transparent theory that accounts for the strongly 



momentum-dependent renormalization experimentally 
suggested in the cuprate superconductors. We focus on 
charge dynamics involved in the low-energy spectra of 
the doped Mott insulators. A key idea is that, near 
the Mott insulator, an electron (a hole) added to an 
empty (a doubly occupied) site costs much smaller en- 
ergy than an electron (a hole) added to singly occupied 
sites and behaves as a component of a band separated 
from the main quasiparticle. Then near the Mott insula- 
tor, this separated excitation contributes to the LUSW in 
the hole-doped (electron-doped) Mott insulators. This is 
in contrast with the weakly correlated regime, in which 
an added electron (hole) constitutes a uniform excita- 
tion irrespective of the added site, because the added 
electron becomes uniformly and spatially extended with 
a good wavenumbcr of the momentum. In our theory, 
such an electron (a hole) added to an empty site (a dou- 
bly occupied site) forms a composite fermionic excita- 
tion (or cofermion), which we call holo- electron (doublo- 
hole). The conventional quasiparticles may be scattered 
by the doubly occupied site (empty site) and are trans- 
formed to the cofcrmions, which introduces finite life time 
into the quasiparticles. This scattering and transforma- 
tion may alternatively be formulated as the hybridization 
between the conventional quasiparticles and the holo- 
electrons and/or doublo- holes. This hybridization nat- 
urally generates a hybridization gap and the resultant 
band splitting causes the Fermi-surface reconstructions, 
namely, topological changes in the Fermi surface. The 
hybridization gap and the topological change account for 
the experimentally suggested pseudogap and Fermi-arc 
formation, respectively, in the cuprates. Our unique pre- 
diction is that the pseudogap is independent of the pre- 
cursor of the Mott gap itself, while the pseudogap has the 
structure of the s-wave type symmetry rather than the 
d-wave symmetry. In addition, the pseudogap identified 
here is in clear distinction from the superconducting gap 
as well. We propose several experimental tests to prove 
or disprove our theoretical consequences. 

A way to understand topological changes in the Fermi 
surface has recently been proposed to originate from 
the emergence of zeros of single particle Green's func- 
tions, which are the points in the (k, w)-spacc satisfy- 
ing Re G(k,uj) = for the momentum k and the fre- 
quency uj. In other words, the single particle self-energy 
T,(k,u>) diverges at the zeros of Green's function. The 
idea of the emergence of the zeros has a root in the 
work by Dzyaloshinskii, who examined a possible ex- 
tension of the Luttinger theorem into Mott insulating 
states^. As is reviewed in Refl29l. this idea is applied 
to the doped Mott insulators to explain the "Fermi arcs" 
observed by ARPES. Recent results of numerical calcula- 
tions also suggest that the zeros emerging in doped Mott 
insulators reconstruct the Fermi surface and changes its 
topolog y 23 ' 25 ' 26 . Such a topological change is indeed 
claimed in a recent ARPES measurement 3 ^. 

The Fermi surface reconstruction itself is not an un- 
conventional phenomenon, if it accompanies a sponta- 



3 



neous symmetry breaking such as an antiferromagnetic 
order. In translational-symmetry broken phases with an 
ordering wavenumbcr Q, electrons with momentum k hy- 
bridize, at least, with electrons at momentum k + Q. 
Then Fermi surface reconstructions are naturally under- 
stood as a consequence of the hybridization gap as we 
will discuss later in detail. It is indeed able to account 
for the formation of Fermi pockets or arc-like Fermi sur- 
faces if a hybridization with other Fermionic excitations 
exist. 

However, in the hole-doped cuprates, in spite of recent 
reports on time reversal symmetry 3 ^ - — and rotational 
symmetry 3 ^ breakings in some of the cuprate supercon- 
ductors, translational symmetry breakings have not been 
universally observed 3 ^. On the other hand, the "Fermi 
arch or the truncated Fermi surface have universally been 
experimentally observed in the underdoped cuprates. To 
have a unified picture of the Mott physics, it is important 
to understand whether the Fermi surface reconstruction 
occurs as a consequence of the hybridization gap gener- 
ated by the "hybridization" with some hidden fermionic 
excitations without assuming a symmetry broken phase. 
In the present theory, the hybridization between the con- 
ventional quasiparticles and the emergent fermionic exci- 
tations called the holo-electron and the doublo-hole nat- 
urally explain such a Fermi-surface reconstruction. 

Now we go into some details regarding characteris- 
tic low-energy spectra of doped Mott insulators. As is 
proposed by Meinders, Eskes and Sawatzky2£, the dop- 
ing dependence of the LUSW formed just above fj,, has 
a marked feature for nearly Mott insulating, correlated 
electron systems. As we have mentioned above, an elec- 
tron added to a singly occupied site inevitably costs the 
on-site Coulomb interaction U and gives an excitation in 
the UHB, whereas electrons added to empty sites exclu- 
sively contribute to the low-energy states just above /x. 
Here, an empty site can accept an electron irrespective of 
its spin state and creates two unoccupied states, namely 
up and down spin states. Therefore the holes cre- 
ate Ah empty sites in the atomic limit, where the kinetic 
energy of electrons t is zero. Then, they create 2N\ l unoc- 
cupied states as a consequence. When the kinetic energy 
of the electrons t becomes nonzero, the number of un- 
occupied states increases because of the pair creation of 
doubly occupied and empty sites induced by the hopping 
of electrons and resultant charge fluctuations. There- 
fore, for hole-doping concentration x = Nh/N s with N s 
being the total number of lattice sites, the LUSW in the 
doped Mott insulator is always larger than 2x. From the 
early stage of the studies on the cuprates, the LUSW 
developing faster than 2x has been observed by opti- 
cal conductivity measurements^. This quick increase of 
the spectral weight larger than 2x requires a picture of 
doped holes very different from the doped semiconduc- 
tors, where the LUSW is trivially equal to x. Our scheme 
presented in this paper naturally explains this unconven- 
tional feature. 

The mechanism of the superconductivity itself is a cen- 



tral open issue of the physics of the cuprate superconduc- 
tors. When the LUSW in the normal state belongs to 
an unconventional phase with emergent excitations, the 
mechanism has to be understood based on this frame- 
work, because the energy scale of the superconductivity 
is even smaller than and governed by the energy scale 
of LUSW. Our theory offers an unconventional channncl 
of the pairing and resultant superconductivity emerg- 
ing from the contribution of composite fermions never 
considered in the literature. We examine an unexplored 
type of quasiparticle pairing arising from the pairing po- 
tential generated by pairing fields of cofermions, holo- 
electrons/doublo-holes and quasiparticles. Although the 
pseudogap formation by the hybridization gap of quasi- 
particles and holo-electrons/doublo- holes is destructive 
to the superconductivity, this unconventional pairing po- 
tential serves to creating superconductors. This dual 
character and two sides of the same coin naturally ac- 
counts for the recent puzzle under debates on the di- 
chotomy and nature of the gap in the anitinodal region 
of the underdoped cuprate o 38 i 39 . 

The organization of this paper is the following: In 
Sec. HH we start from the Kotliar-Ruckcnstcin's mean- 
field theory and review previous extensions for correlated 
metals for the self-contained description. In Sec. IIIII we 
take into account the charge dynamics, which plays an 
important role in formation of the LUSW and explain 
how emergent excitations, namely, the holo-electrons and 
doublo-holcs emerge. The hybridization between quasi- 
particles and these composite fermionic excitations nat- 
urally causes the Fermi-surface reconstruction, namely, 
the topological change in the Fermi surface. The pseu- 
dogap phenomena observed in cuprate superconductors 
also emerge because of this hybridization. In Sec. IIIID . 
we propose a novel pairing mechanism evolved from the 
cofermions as the mechanism of the high temperature 
superconductivity. Our results for single particle spec- 
tra, amplitudes of the pseudogap, Fermi-surface topol- 
ogy, the specific heat coefficient, and the density of states 
are presented and compared with experimental results in 
Sec. HV] We also estimate the superconducting gap am- 
plitude and quantitative aspects of the superconducting 
mechanism. Sec. [V] is devoted to discussions and sum- 
mary. 



II. PREVIOUS THEORIES 

A. Hubbard model 

The Hubbard mode l 15 ' 17 ? 40 defined by the following 
Hamiltonian, 

is a canonical model, which describes the competition 
between the kinetic energy tij and the on-site Coulomb 
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interaction U, where cJ CT (c icr ) is a creation (annihilation) 
operator of cr-spin electron on i-th site and h ia = c ia c ia 
is a number operator. Hereafter we focus on the Hubbard 
model on the square lattice as a model for the cuprates. 
In this chapter, the nearest-neighbor hopping and the 
next-nearest-neighbor hopping are set as iy = — t and 
tij = +<', and further neighbor hoppings are ignored. 

The solution of the Hubbard model on two-dimensional 
lattices remains an open problem. To get an in- 
sight into the nature of the correlated metallic phase, 
there exist a variety of numerical methods, which give 
accurate results for finite size clusters such as exact 
diagonalization 4 ^ and quantum Monte Carlo^. Cluster 
extensions of the dynamical mean-field theor y 43 i 44 , im- 
proved variational Monte Carlo method s 45 ' 46 , Gaussian 
Monte Carlo metho d 47 ' 48 and the path integral renor- 
malization grouper— are also available in the literature. 
However, they all have some limitations. For example, 
limitations on cluster size are severe in the exact diago- 
nalization and the quantum Monte Carlo, while resolu- 
tions in the momentum space is severely limited in the 
cluster extension of the dynamical mean-field theory. 

Without relying on accurate numerical methods, there 
exist other ways to extract essential physics from analyt- 
ical or conceptually correct limits. The Landau's Fermi- 
liquid picture 5 ^ offers such an example. Although we 
have no exact solutions in the thermodynamic limit, low- 
energy single-particle spectra of the Hubbard model are 
believed to behave following the Landau's Fermi-liquid 
picture for less correlated systems except for ID systems. 
The original proposal for the Mott insulating states itself 
is not originally based on the numerical results, but pro- 
posed by Peierls and Mott from a gedanken experiment 
on metallic crystalline hydrogen-like atoms^r— . Such 
phenomcnological theories based on physical intuitions 
offer insights into difficult issues in a wide range of con- 
densed matter physics. In this paper, we try to make a 
step towards constructing such a physically transparent 
theory that accounts for the unconventional properties of 
doped Mott insulators. 



We start with the local Hilbert space of the Hubbard 
model, which is expanded by a set of the Fock states; 
the empty state (holon) |0), the singly occupied states 
| t), | I), and the doubly occupied state (doublon) | f! 
) . Corresponding to each Fock state, one slave boson 
is introduced as e for |0), p a for |cr), and d for | fi 
) . In addition to these bosons b (b = e,p a or d), a 
fcrmion operator f a is introduced to stand for the cr-spin 
quasiparticle. The correspondence relation to local basis 
for the lattice wave functions is given as 



\0)i 
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t ^|0)i = /M|vac; 



(2) 



(3) 



(4) 



(5) 



where |vac)f^ B ' 1 is a vacuum of the i-th site for fermionic 
(bosonic) degrees of freedom. Equations ©-(O represent 
the mapping between wave functions written by the orig- 
inal electrons and the fermions / combined with bosons 
b. This mapping is derived when the electron operators 
are replaced with composite ones as 



(6) 



We should note that the replacement given by Eq.© is 
not a unique one. It is known that operators equivalent 
to the right hand side of Eq. © can be given as 



ct = z- ft 



where z,„ is define d 56 ' 57 as 



9ia \p\ak + djfc) gfj, 



(7) 



(8) 



B. Kotliar-Ruckenstein formalism 

One of the simplest picture to describe correlated met- 
als and metal-insulator transitions at half filling (Mott 
transitions) along the line of the original idea by Mott is 
the Brinkman-Rice scenario^ 6 - based on the Gutzwillcr 
approximation^ 7 -. As our starting point, we employ 
the Kotliar-Ruckenstcin's (KR's) slave-boson formalism, 
which gives the same results as the Gutzwillcr approxi- 
mation. In this subsection, we briefly outline the KR's 
slave-boson formalism for the Hubbard modelS 6 -, to make 
the present paper self-contained. This slave-boson for- 
malism gives a starting point for the mean-field descrip- 
tion of strongly correlated electron systems by replacing 
original electrons with four kinds of bosons and one kind 
of spinful fermion. 



€ } = (l-4fe-4eV) P \ (9) 
€ ] = ^-v\j ia -d\d,) P2 . (10) 

The operators g£ and act as identities when these 

are operated to (j>\ tT e i + d\p— ^ , for any powers pi and 

P2- This ambiguity of the correspondence has been uti- 
lized before^ 6 -. 

In the expanded Hilbert space, the Hubbard Hamilto- 
nian Eq.([T]) is rewritten as 

H = J2 t^JUiJl + uY,d\d z . (11) 

The on-site Coulomb interaction is replaced by a "chem- 
ical potential" for doublons d\. Then the correlation 
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among electrons is now contained as hopping process of 
quasiparticles fj a disturbed by the associated motion of 
slave bosons. The quasiparticle hopping is accompanied 
with four kinds of bosonic motions generated by z ia Zj a , 
namely, physical processes of hopping of holons ej and 
doublons eft, pair creations and annihilations of holons 
and doublons. 

When we employ the path integral description of the 
system by making use of the coherent states for bosons 
and fermions, we need to introduce a set of constraints 
to eliminate unphysical states in the expanded Hilbert 
space, which arise when we introduce slave bosons to de- 
scribe the local Fock states. First, only one boson should 
occupy each local state. There are only four local phys- 
ical Fock states, and these four states are exhausted by 
four different kinds of slave bosons. Therefore, we need 
the first constraint 



E 2 ^ 



1. 



(12) 



In the path integral form, the partition function for the 
Hubbard model is given by 



Z = Jv[clc a ] 



,~S\cl,c„ 



(14) 



where the action is 



S [4 , c CT ] =Y,f dT 4to [(fir - M) Sa + t tj ] c ja {r) 
+ UJ2 / dr cj^c^cl^c^r). (15) 



We use the same notation for both of the operators and 
the corresponding grassmann fields for simplicity. 

On the other hand, the partition function in the slave- 
boson formalism is 



Second, the number operator of the cr-spin quasiparticle 
is necessarily given as 



f T f 

J 1(7 J 11 



PiaPia 



did,. 



(13) 



Z= IV fl /J V \b\b\ V [A] e- s . (16) 



The action is rewritten as 



S = J drflir) [(d T -n + \%) 



kj+QMUj fM+ dr UdK^d^r) 



ir £ ej^fire^r) + $>l(r)d T p i(T (T) + cI1(t)cU(t) 



(17) 



where Si = {e\,p\^,p\^,d\), A, = (A,- , X^r, X\f), and These saddle point values are determined sclf- 
l ( \ _ 2 ( \s1 ( \ Th LaeranEe multipliers A (1) consistently through minimizing the free energy / 

given in Eq.® below. Then the action for the 
fcrmionic degrees of freedom contains only quadratic 
terms of fermionic fields j\ G after the slave bosons are 
replaced with c- numbers eo, Po a , and do as 



(2) (2) 

X\J , and X\ , are introduced to reinforce the set of con 
straints Eqs. $TZ§ and (jig]). 



C. Mean-field theory for KR formalism 

The mean- field approximation for Eq. (|16[) corresponds 
to replacing bosonic fields Bi and A^ with the homoge- 
neous saddle point values for them as 

ej, e i ->• e , p\ a , p la -^P 0a , d\, d t ->do, 



^E fl^ ^ - >*) + Co**y] /,>(r).(18) 



We can easily integrate out the remaining fcrmionic de- 
grees of freedom and obtain the mean-field free energy 
for homogeneous phases as 

/ = -IL^l n [l + e -' 3 (^-^ 2) )" 



Ud n 
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(19) 



where N s is the number of sites, efc is the Fourier trans- 
formation of tij and the mean-field quasiparticle renor- 
malization is given as 



Co* = [gia {Po* e O + d p ow ) g 2 a] ' 

/1 —2 — 2\P1 
9ltr = (l-POa - e 0j J 

1 -P0o 



(20) 
(21) 

(22) 



It is known that with this drastic mean-field approxi- 
mation, the non-interacting limit is correctly reproduced 
when we set p\ and pi defined in Eqs.© and (|TU)) as 
Pi = P2 = —1/2. For [7 = 0, the mean-field val- 
ues of the density of bosons are given as, d = (n/2) 2 , 

~Eq = 1 — n + d , and = p^ = n/2 — d . Then the 
mean-field renormalization factor turns out to be 1, cor- 
rectly. Hereafter, the values of p\ and pi are fixed as 
Pi =P2 = -1/2. 

The saddle point values for bosonic fields, eo, p 0fT , 
and do, can easily be examined in some well-defined lim- 
its. In the strong coupling limit U/\t\ ^> 1, the density 

—2 

of doublon is suppressed, d ~ 0, at any doping level. 
Then the spin density is roughly proportional to the spin- 
dependent electron density, p% a ~ n„ (please recall the 

—2 — 2 \ 

mean-field constraint, n* = p Qa . + d ). In the hole-doping 
case, the density of empty site eg is given by the doping 
concentration, e§ ~ x. 

In this limit, the doping dependences of (o* and LUSW 
are given as follows: By using the fact n a — n/2 = 
(l—x)/2 in the paramagnetic phase, the mean-field renor- 
malization factor £o is given as 



Co* — Co — 



Po ( e o 



d y 



l-x 
2 



2 



2x 



l+x 
2 



1 



(23) 



Multiplying the number density of unoccupied states, 1 + 
x with Coj we obtain the LUSW as 2x. This doping- 
dependence of the LUSW is the same as that of the exact 
solution in the strong coupling limit. 



D. Previous studies on charge fluctuations 

1. Formation of upper and lower Hubbard bands 

The mean-field theory only accounts for the coher- 
ent quasiparticle excitations in the correlated electron 
systems. Of course, momentum-dependent renormaliza- 
tions do not appear. For overall description of the energy 



spectrum including incoherent Hubbard bands, dynam- 
ics of the charge bosons is known to be essential when 
one wishes to improve the slave boson mean-field the- 
ory. For example, Castellani et al. have claimed that the 
Gaussian fluctuations of charge bosons around the saddle 
point solution can reproduce the structure of incoherent 
Hubbard bands^. To take into account the fluctuations 
of bosonic fields around the saddle point solutions, the 
Bogoliubov prescription^ is used, in which the boson 
operators arc divided into condensate components and 
fluctuating components out of condensation as 



et = e + e|, e i 



eo 



Pi* = PO* + Pi , Pi* = PO* + Pi , 

d\ = d + d] , di = do + ds ■ 



2. Bosonic propagators 



(24) 
(25) 
(26) 



Propagators of the Gaussian fluctuations of the slave 
bosons are given by the action Sb with quadratic terms 

including only the bosonic fields B . B as 



j(0) c (ic) c (is) _,_ g(lcs) 



(27) 
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f dr ^{el^^ + Af^e^r) 
Jo 

d}(T)[d T + U + X^-J2^J]di(r)}, (28) 

cr 

f drY J ^i{r)Lt j ^{r), 
J ° ij 

[ P drJ^^irW^Mrf 
Jo 

a 

dTj2[Pi(T)L1» a cl> jt T (T) 



(29) 



(30) 



+ 4j(r)I^J^{r) + 6 l3 KZ 

x {^(t) ■ ct>l(T) + 4> ta (T) ■ (3l(T)}} , (31) 

where we use vector notations, j3 i = (e i ,d} i ), <p i = 
(p icr ,pt-), and coefficients defined as 



-c . ( P 2 o* +pI* 1Po*Po* 

—2 



e + d 2e d 
2e d e + d 



(32) 
(33) 
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(34) 
(35) 



Here the hopping parameter for bosons is given as <7y = 
^j'fl'icrfl'Icr (fiafj<r/i wncrc static correlation functions for 
quasiparticles (^f^fjj^ arc introduced. The average (• • •) 
is defined as 



(X) = 





j a i J a 


X e -So-s B 




f<7i fa 





(36) 



where Sq + Sb is the approximate action used in this 



subsection. The term \fl a fj a j represents the mean 

field that stands for kinetic motions of quasiparticles 
surrounding bosons. This mean field seems to be 
self-consistently determined through Eq.([3l)|). However, 
by using the approximate action, Sq + Sb, we obtain 
ft f. 

J ujJ ja 



ft f ■ ) = 

J %a J ja j 





fa J fa 


ft f. e -So 
J zaJ ja^ 




fai fa 


e -S 



(37) 



because Sb does not contain f ia and/or fj a . This mean 
field (firfj,/) is obtained through decoupling the "inter- 
action" term z la f}J ja z] a as 

%iafiafja Z ja — ^ia^ja (fiafja/ + (^ia^ja/ fiaf 



^■>a^ja ) {fiafja 



(38) 



by neglecting fluctuations 



3° 



V f - < V f 

JiaJja XJzaJja 



Averages such as (f^fjj) and (m^z^JS are taken by 
using the resultant action. However, we should note that, 
in the previous studie a 57 ' 58 , the average (z ia z^J\ were 
treated as £o and the contribution from static correlation 
functions of bosons such as ^dje^\ was dropped. We 

also note that, in the mean- field level, the dispersion of 
spin bosons for the paramagnetic phase vanishes as x — > 
0, because of simultaneously vanishing condensation of 
charge bosons: eo, do — > in Eqs. (|30| and (|3"Tj). 

In the close proximity to the Mott insulating states, 
where So, d <C 1, propagators for the Gaussian fluc- 
tuations of charge bosons f3 i = (/3j,(3f) = (e i7 3\) are 
approximately determined by the action S^ + Sg 1 ^ as 



0V 



6t \ _ z+(Q) 



Q 



Z a _ b {Q) 



iu) m -\K Q \ iu m + \\ Q y 



(39) 



where the coefficients Z±(Q) are given by 



z±(Q) z?(Q) 

Zf{Q) Zf{Q) 



1 



SX + SU/2 ( l 
2a Q \ 1 



' —1 J 2a Q \ 1 1 

Parameters used in the above equations are given as 



(40) 



A (1) - A (2) , 



(41) 



^ 1) + ^-vIWq) -Pokl 2 4' ( 42 ) 




su 
T 



(43) 
(44) 



where SU = U-2X^ gives the amplitude of the Hubbard 
gap at half-filling n = 1, and 



1 



Raimondi and Castellani^ introduced the following 
approximate form of the single particle Green's function 
as 

GijAr) = -(Tc iCT (r)ct CT (0)) 

= -(Ti(r)z 3ff (0)/ lff (r)4(0)) 

* - (4(t)V(0)> (TfMfUo)) ,(45) 



where vertex corrections are dropped. A further ap- 
proximation was introduced: By focusing only on the 
charge dynamics, contributions from the fluctuations of 
spin bosons were neglected as 



(t4(t)^(o)) ~ Coa+pfolglo 



X (Tb ia (T)bl(0) 



(46) 



where b irT (r) = e| (T (r) + d ia (r). As a result, the approxi- 
mate Green's function is given as 

g ij<7 ( T ) ~ -to*(Tf ia (T)fUa)) 

- pl 9 la9la{T\a{r)b]a^)) 

x (TfMfUo)) . (47) 



The first term of the right hand side of Eq. (|4"T|) gives the 
coherent band and the second term gives the incoherent 
Hubbard bands. 
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III. COFERMION THEORY 



A. Perspective 



eifc 



A 2 

AFM 



£2fc + M + V (\ k + A 



2 

AFM 



(49) 



In the above sections, we reviewed the mean-field KR 
theory and how fluctuations of charge bosons induce the 
incoherent bands. These previous theories have draw- 
backs, in spite of an advantage in the simplicity. For ex- 
ample, the momentum-independent quasiparticle renor- 
malization in the previous theories cannot account for 
the "Fermi arc" observed in ARPES measurements of the 
cuprate superconductors^. Since significant momentum- 
dependent quasiparticle renormalizations have been cap- 
tured in numerical works through the cluster type exten- 
sion of the dynamical mean-field theory**—, it is desir- 
able to construct a theory that is physically transparent 
and can examine the observed singular self-energy, which 
results in, for examples, the pseudogap and the Fermi 
arc. We assume that the experimentally observed arc- like 
Fermi surface, in the hole-underdoped cuprates, is a con- 
sequence of Fermi-surface reconstruction caused by the 
divergence of the quasiparticle self-energy or emergence 
of zeros of the Green's function. The divergence means 
the breakdown of the perturbation theory and hence the 
breakdown of the Fermi liquid theory as well. Here we 
extend the previous theories to understand this uncon- 
ventional feature. Our goal is to acquire a simple frame- 
work and an intuitive understanding. 

The Fermi-surface reconstruction itself is not an un- 
conventional phenomenon as we have discussed in Sec. U 
in the example of the antiferromagnetic order in the ordi- 
nary Slater's mean-field description. In an antiferromag- 
netic metal on square lattices with the ordering vector 
Qi) = (7r,7r), the f-spin electrons with momentum k hy- 
bridize with the ,J,-spin electrons with momentum k + Qo, 
as the term AAFM[cL.c fc+ g o ^ + H.c], in the presence of 
the mean field Aafm- As a result, the self-energy for the 
t-spin electron diverges at the momentum k, where the 
^-spin electron with momentum k + Qq has a pole. 

When the bare band dispersion of ct before the hy- 
bridization is switched on is given by 

£k = eifc + £2fe — Mi 
eik ~ —2t(cosk x + cosk y ), 
£2k = 4t' cos k x cos k v , 



the Green's function of the electron <^ ka that hybridizes 

k+Qoa 

1 



with the electron ct , ^ — is obtained as 



G(k,u) = 



A 2 

^AFM 

1 - ffc+Qo 



(48) 



2 V*ffc + A AFM 



The Green's function written in Eq. (|48l) shows that 
G(k, w) vanishes at w - Cfe+Qo = 0. I n other words, 
G(k,ui) has zero at uj — Cfc+Qo = 0- This emergence 
of zeros is nothing but the divergence of the self-energy 



of the electron c ka , E(fc,w) 



a Ifm/( w 



Cfe+Qo) at 

w — £fc+Q = 0, as is seen in Eq. (|48| . As a result, 
the zero surface defined by uj = £,k+Q a splits the pole 
surface defined by uj = |fc int o the two pole surfaces 
u = e 2 fc - /i ± \/ej k + A afm> as is seen in Eq.(|49]). 
Then the Fermi surface is disconnected into pockets at 
the gap edge. As an example, a zero surface, recon- 
structed band dispersion, and Fermi surface are depicted 
in FigH for t' = 0.25*, \i = -0.75*, and A A fm- This 
has to some extent a qualitative similarity to what is ob- 
served in the cuprates. Another example is the case of the 
BCS superconductivity where the quasiparticle cj^ has a 
"particle-particle hybridization" proportional to cL.c_ fe ^ 
with cL fc4 .— ■ 



tii=0 (n,n) 




UJ 



£2k + i" - x/ef 



lk 



"a 2 " — 

AFM 



FIG. 2: (color online) Left panel shows band dispersions and 
zero surface along lines running from (n,n) to (0,0), from 
(0,0) to (ir,0), and from (tt,0) to (0,0). Right panel shows 
bare and reconstructed Fermi surface, and zero surface at 
uj = 0. Thin (black) solid curve stands for bare band dis- 
persion u — £fc , and thin dashed (black) curve stands for zero 
surface u> = £fc+Q - Bold (red and blue) solid curves stand for 
reconstructed bands. 

However, as is already remarked in Sec. U symme- 
try breakings such as the antiferromagnetic orders have 
not been universally observed in the hole-underdoped 
cuprates where the "Fermi arc" or the truncated Fermi 
surface has been experimentally suggested. Therefore, 
we need a mechanism for the emergence of zeros with- 
out any symmetry breakings. In this context, the trun- 
cated Fermi surface scenario reviewed in Refl29l gives an 
interesting example inspired by a result of the renormal- 
ization group methods, which attributes the emergence 
of zeros to Umklapp scatterings without the symmetry 
breakings^. The pseudogap in this framework is the 
precursor of the Hubbard gap. We employ an alterna- 
tive framework, where a zero surface with a gap for- 
mation naturally emerges by a hybridization with some 
additional fcrmionic excitations, as is discussed in the 
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above example of the antifcrromagnetic ordered phase, 
although our additional fermions arc different from the 
quasiparticlc at k + Qo in the antifcrromagnetic case. 
Since our pseudogap (the hybridization gap) will turn out 
to be different from the remnant of the Hubbard gap, our 
framework will yield results qualitatively different from 
the scenario by the Umklapp scattering as we see in this 
paper. 

Is there such fermionic excitations in proximity to Mott 
insulators in the absence of symmetry breakings? Or how 
do such fermionic excitations emerge? A hint for the exis- 
tence of such additional fermionic excitations comes from 
examinations of the LUSW illustrated in Fig. [TJ in the 
doped Mott insulators. Hereafter we restrict our discus- 
sion to the hole-doped Mott insulator as in Fig. [1] with 
the hole-underdoped cuprates in mind. The LUSW is de- 
fined as the spectral weight above [i within the coherent 
band near the top of the LHB. 

We first recall the origin of the LUSW discussed in 
the literatures^. First we begin with the atomic limit, 
where t = and U ^ 0. Then the LUSW of the hole- 
doped Mott insulator consists only of states created by 
adding an cr-spin electron to an empty site, namely, a 
holon site, to avoid creating a doubly occupied site, dou- 
blon, and to avoid the cost of the on-site Coulomb repul- 
sion U. An electron added to an empty site is confined 
tightly, in this limit. If the electron escapes from the 
holon sites, it inevitably creates a doublon and costs U. 
In other words, an electron created at a holon site can 
contribute to the LUSW, although this electron does not 
propagate coherently. 

In the KR formalism, creation (annihilation) operators 
for the electron at the holon site are given by composite 
fermion operators e i fj a (/j CT eJ)- This is just a cr-spin 
electron in the original Hubbard model and is definitely 
fermionic. 

When t/U becomes nonzero, we have tightly bound 
doublons and holons even in the Mott insulating phase. 
Therefore, holons in the hole-doped systems consist of 
both of the doped holons and the preexisting holons al- 
ready present in the insulators. In the hole doped sys- 
tems, these two types of holons should not be distin- 
guished and should constitute the same object. Then an 
electron added to these holons should constitute LUSW 
and may constitute a novel composite particle distin- 
guished from the original electron and offers a hint for 
the additional fermionic excitations which brings about 
zeros to the quasiparticlc Green's functions, if this com- 
posite particle hybridizes with an original quasiparticlc. 

Originally the binding energy of a doublon and a holon 
in the Mott insulator is the order of U . However, for finite 
t and finite doping x, indeed, quantum fluctuations in- 
duced by coherent carriers dramatically weaken bindings 
between doublon-holon pairs. Then an electron added to 
this holon bound to a doublon requires only small en- 
ergy and merges into the excitation of an added electron 
to a doped holon site. Actually, remnants of doublon- 
holon pairs, namely, weakly-bound doublon-holon pairs 



are known to play an important role in correlated elec- 
tron systems, especially in the context of variational wave 
function thcorie o 63 ! 64 . 

A part of the tightly-bound doublon-holon pair, eof} a , 
has already been taken into account in the mean-field 
level in the previous theories (please note e i f} (J = eo// CT + 
^iflcr)- However, it represents nothing but the renor- 
malized quasiparticlc, which propagates in homogeneous 
mean fields and has nothing to do with the above com- 
posite particle. To substantiate our picture, we need 
to take into account the composite operator including 
bosonic fluctuations, e i fj !T . If we treat such a compos- 
ite operator e i f\ a as a fermion, overlaps between tightly- 
bound doublon-holon pairs and quasiparticle states cause 
a hybridization of fermions between two types, where one 
is the quasiparticle f\ a and the other is the "composite 
fermion" or "cofermion" e i p ia . From such a hybridiza- 
tion between the cofermion and quasiparticle, weakly- 
bound pairs discussed above will naturally be born as a 
result. 

If such a hybridization between the quasiparticle and 
the composite fermion really exists, this hybridization 
would contribute to the self-energy of quasiparticles. De- 
pending on the dynamics of the composite fermions, such 
self-energy would have a strong momentum-dependence 
and possibly have poles near the zero energy In such 
cases, the hybridization between the quasiparticles and 
the composite particles causes a hybridization gap near 
the Fermi level. This gap is expected to account for the 
pseudogap phenomena. As is discussed below, we con- 
struct an action containing both of the quasiparticles and 
the composite fermions. 

In the following sections, we concretely give our theo- 
retical treatment by introducing the composite fermions 
discussed above. First, we list up shortcomings of the 
previous theories and present our guiding principle to 
overcome the previous KR's mean-field theory. 

As is mentioned in the above section, the KR formal- 
ism gives exact results if we thoroughly take into account 
the fluctuations of slave bosons and Lagrange multipli- 
ers introduced to keep the constraint Eqs. (fT2]) and (|T3| ( 
Hereafter we call this constraint the local conservation ) . 
However, the fluctuations of slave bosons are neglected 
in the original KR's mean-field theory. In addition, the 
constraint is also treated just on average. The original 
constraint is for the operators or the dynamical bosonic 
fields. In the mean-field theory, the constraint is replaced 
with averaged density of boson condensations. When we 
include the fluctuations of slave bosons, we should take 
care of the constraint or the local conservation of the 
densities of the slave bosons. In the previous theories, 
by including the fluctuations of charge bosons, the inco- 
herent Hubbard bands are reproduced. However, there 
exists a difficulty in conserving spectral weight, because 
the fluctuations violate the local conservation of the bo- 
son density^. 

When we impose the local constraints more strictly for 
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fluctuating bosons beyond the mean-field level, it turns 
out that the term 

Qjl = 9la9la (pL^i + APia) (^jPja + P^j ) ( 50 ) 

represented by the diagram in FigJHJa) is dominating 
among all the possible diagrams for the kinetic term in 
the action Eq. (fTT]) . namely 

/ dr £ ft {r)tdija (r)f) a (r). (51) 
Jo 



Here we employ g\ a = (1 - pg- - eg) 1 and g\ a = (1 
Poct — ^o) _1 by following ReflBH 



(b-2 



(C-2) 




(h) 



FIG. 3: (a)-(e): Diagrams representing terms in fJ a Cijfj a - 
(f)-(h): Examples of various diagrams on the one- loop level 
generated from the coupling (polygons) and derived from the 
terms (a)-(e). Solid lines with arrows represent propaga- 
tors of the quasiparticles. Wavy lines stand for the charge 
bosons, and bold solid lines are the spin bosons. Condensa- 
tions of bosons are represented by lines terminated at crosses. 
Coupling constant gi a g2a-tij is represented by open polygons. 
Here, we do not distinguish holons and doublons. Spins are 
also not distinguished in the diagram, for simplicity. 

To elucidate why we retain Qjli we classify the dia- 
grams illustrated in Fig |3] into four types, categorized by 



the time dependence (or frequency dependence) of quasi- 
particles and the local conservation of the boson densi- 
ties: 

(T-l) Diagrams containing external propagators of 
quasiparticles, in addition to bosonic propagators 
violating the local conservation of boson densities 
(Figures [3]3 f3j; J3ji , and [3ji) . Here the violation 
means that before and after the interactions (rep- 
resented by hexagons), the number of bosons ex- 
pressed by external boson propagators is not the 
same. 

(T-2) Diagrams that contain time dependence of quasi- 
particles, but that do not violate the local conser- 
vation (Figures 13^. and [St:). 

(T-3) Diagrams that do not contain time dependence of 
quasiparticles but do violate the local conservation 

(FigEF). 

(T-4) Diagrams that neither include time dependence 
of quasiparticles nor violate the local conservation 
(FigEfe). 

Here we present our guiding principle to take account 
of boson fluctuations: We exclude (T-l) because it vio- 
lates the local conservation when the quasiparticles dy- 
namically fluctuate. On the other hand we retain dia- 
grams belonging to the categories (T-2), (T-3), and (T- 
4). The reason to retain these diagram is as follows. 
The diagrams in the category (T-2) do not violate the 
local conservation, when bosons fluctuate. Therefore, we 
take the diagrams in this category into account. On the 
other hand, the diagrams in the category (T-3) do vio- 
late the local conservation. However, in these diagrams, 
quasiparticles enter as time averaged Green's functions. 
Therefore, quasiparticles feel the time averaged bosonic 
motions. The real violation of the local conservation oc- 
curs only when a dynamical quasiparticle process is in- 
duced by fluctuating boson hoppings. On the contrary, 
the real violation does not occur when the quasiparticles 
emerge as the time averaged quantities as in the case of 
(T-3). This is the reason to retain the diagrams in the 
category (T-3). Since (T-4) does not violate local con- 
servation, we retain it. 

For the slave-particle formalism of correlated fermion 
systems, it is wellknown that fluctuations of gauge fields 
play an important role on reinforcing the local constraint 
imposed on slave particles^. It was pointed out by Joli- 
coeur and Le Guillou that the Kotliar-Ruckcnstcin for- 
malism has the U(l)xU(l)xU(l) gauge symmetr y££. It 
comes from the phase symmetry of the slave bosonic par- 
ticles, namely, e,, p ia , and d i . 

In our theory, we will treat fluctuations of such 
phases together with fluctuations of the amplitude of 
the condensation fraction of these slave particles, by 
using the Bogoliubov prescription. Therefore, the 
phase fluctuations are taken into account, although the 
U(l)xU(l)xlI(l) gauge structure is not strictly con- 
served. 
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B. Stratonovich-Hubbard trasformation 

We introduce Grassmannian valuables (or fermionic 
fields) T ia = (i/J ia , Xia) T that stands for the cofermions 
as are discussed in conceptually in Sec. IIIIIA . by using 
the following identity 



(a) 



]drldr t<J e A = dct 



where matrices T a are defined as 



{T a )ij — .9l<j32a^ 



(52) 



(53) 



A = f" drJ2[(rUr)-Cl(i 

ija 



and 



x {Ta)ij ^T jo .(t) — C J(T (r) 
Here wc use vector notations as 



Ci a — {e i ,d}^j fj al C ia — f ic7 (e\,d i 



(54) 



(55) 



The identity Eg. (f5"2"]) gives the transformation for 
a coupling term of the quasiparticles and fluctuating 
bosons depicted in FigJ3^, 



(56) 



cxp 
where 



drC a 



J ia 

dot [T t Tj 

a = £■ al + £■ a2i 

C' al = Y^r\ a {T){T a )i^ j(T {r), 

iju 

C'a2 = -^{cL(r)(f ff ) i: -f 3V (r) 

ija 



(57) 



(58) 



(59) 



(60) 



These transformed Lagrangian C' a i (FigHh) and C' a 2 
(FigHb) lead to the cofermions' self-energy and the hy- 
bridization between the quasiparticles and cofermions, 
respectively after integrating out the fluctuating bosonic 
degrees of freedom. It will be discussed below, by using 
a set of the Dyson equations. 



(b) 



FIG. 4: Diagrams for transformed Lagrangians defined in 
Eqs. (|58l) - (|60[) . Solid and dashed lines with arrows represent 
propagators of the quasiparticles and cofermions, respectively. 
Wavy lines stand for the charge bosons, and bold solid lines 
are the spin bosons, (a) The diagram represents the La- 
grangian C' ai fEq.(|59|lV (b) The diagrams stand for terms in 
the Lagrangian £' a 2 (Eq.([5DJ)). 



C. Prescription for self-consistent procedure and 
Green's functions 

Here we construct approximated Green's functions for 
the Gaussian fluctuations of the bosons, quasiparticles, 
and cofermions by using a set of Dyson equations as 
is depicted in FigJS] Wavy lines and bold grey lines 
stand for the Green's functions of the charge bosons 



TPf(T)pf (0)) , and the spin bosons 

Abt 



A ab {r, t) 



C ab (r, t) = - (T^fiT)^ '(0)J, respectively, where a, b = 

1,2, r = i-j, = and = 

(Pia'Pia)- Bold lines with arrows represent the quasi- 
particles propagator Q„ - (r, r). On the other hand, Thin 
wavy lines and thin lines represent bare propagators of 
the charge bosons AQ b (r, t), the spin bosons C^ b (r, r), 

respectively, determined by C^' , in which self-energy ef- 
fects arc not taken into account. Thin lines with arrows 

( f) 

stand for bare propagators of the quasiparticles Q 0rJ [r, r) 
determined by 



C 



fl(r)[(d T - m)^- + (o*U 3 ]f 3a (T), (61) 



where ( 0rJ = g\„giAPoa^ + d oP a) 2 - The Lagrangian £ 
is obtained by decoupling the fluctuating bosons from 
the KR's action (see Eq.(|17|)). Bold and thin dashed 
lines stand for the cofermions' propagators T ab and bare 
propagators = <5 Q ,b/e (e — > 0), respectively. 

In the set of Dyson equations (FigJ5J), wc neglect the 
coupling between charge and spin bosons described by 

propagators such as (p\ a ei \ , at the Gaussian level, since 

these coupling terms are higher order contributions. Be- 
low we explain that the coupling gives higher order con- 
tributions with respect to the hole-doping rate x, in prox- 
imity to Mott insulating states: Since operators including 
both charge and spin such as p\ IJ e i do not conserve the 

electric charge, propagators such as (pl^i ^ should van- 
ish in the Mott insulating phase, where the charge can 
not fluctuate. Therefore, the charge and spin excitations 
are well separated in the Mott insulating phase. 

When hole carriers are doped, jjf-e^ can have a nonzero 
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FIG. 5: Diagrams for Dyson equations. Solid and dashed 
lines with arrows represent propagators of the quasiparticles 
and cofermions, respectively. Other notations are the same as 
FigEl Bold (thin) lines are for renormalized (bare) propaga- 
tors. Filled (blue) circles are amplitudes of the hybridization 
between quasiparticles and cofermions. 



expectation value, at most, scaled by the condensate frac- 
tion of holons eo which gives a rough estimate of the am- 
plitude of charge fluctuations. From a relation e 2 , oc x 



held in the KR theory for the hole-doped case, we obtain 
Jx~. Futhermore, there is an additional con- 



strain! for the coupling terms such as ^pJo-e^ ) : they do 
not appear alone in calculations of physical quantities. 
To conserve charge and spin on average, (pl (T e i \ appears 

with (e\p i(J \ in pair, for example. Therefore, the contri- 



bution of the coupling between charge and spin bosons to 
physical quantities is scaled by (^/x) 2 . It concludes that 
the coupling between the charge and spin bosons gives 
contributions as a higher order in terms of x in physical 
quantities. 

By solving the set of Dyson equations, we obtain the 
propagators for the quasiparticles and cofermions. Here 
the bosonic degrees of freedom are taken into account in a 
self-consistent fashion, through the cofermion self-energy 
(r, r) , and the amplitude Ajj of hybridization be- 
tween the quasiparticles and cofermions, each of which 
we detail below. 

The Lagrangian for the cofermions is given by 



4f =-x;*L(t) 



(r,r) Tj CT (r), 



(62) 



where Y 



(Vkt'XL) i s a vec t° r notation for the 



cofermions, as is defined in the main article, and r = i—j. 
The cofermion self-energy S^ cf ^ (r, r) is a 2 x 2 symmetric 
matrix, 



S (cf) 



n 1 K 



(63) 



The details for the self-energy matrix are given in Ap- 
pendixO On the other hand, the hybridization between 
the quasiparticles and cofermions is described by 

£ hyb = £[f L(r)A^4(r) + &(T)A£f ja (r)}, (64) 



where Ay = (A^ , A\f ) . For details for Ay, see Ap- 
pendix [Bj 

As a result, the effective Lagrangian for the quasipar- 
ticles and cofermions, C c g, is given as 



•£ e ff = £o + -Ccf + £ 



hyb- 



(65) 



When the charge gap is relatively small, E^ ~ E;; and 



A\j ~ A\J ' hold approximately. When the charge gap 
collapses, E* 1 = E 22 and = A[f hold exactly In 

our results, we employ approximate relations E = E* 1 ~ 
E 22 and An = A^ ~ A$. Then a cofermion mode 

vPka + XfcoO/v^ hybridizes with quasiparticles through 
the amplitude A(fc), which is depicted in FigJ5]as closed 
(blue) circles, where k is a momentum. The inverse of 
cofermion propagator (namely, the cofermion self-energy) 

for (VV + Xka)/ is given as 



(x) 



1 



[E a (fc, is n ) + E^fc, ie n )] = lk ie n - a k + 0(4),(66) 



where e n is a fermionic Matsubara frequency. 

Then, the Green's function for the quasiparticles is 
given as 



G ( J\k,ie n ^uj + i5) = G y J>(k,uj + i5) 



A(fc) 2 



7fc(w + iS) - a k 



- (67) 



where e k is the Fourier transformation of t\j , and y, is the 
chemical potential. Here we note that the weights of the 
two quasiparticlc bands split by the zero surface defined 
by uj = a k /lk are not the same in our theory. 

In our calculations, we define the doping rate x by 
using the quasiparticle Green's function as 



T V- 

1 — x = lim — > 

T^0+ N s ^ 

k,i€n ,C 



g { J } (k,ie n ), 



(68) 



where T stands for temperature and N s is the number of 
sites. 
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The Green's function for the electrons, instead of the 
quasiparticles, is given as 

GijAr) = -(t^(t)c] ct (0)} 

~ -(t4(t)% ct (0)) 

x (T/, CT (r)4(0)) 

= (r4(r)V(0))^(r), (69) 

where the bosonic and fermionic degrees of freedom are 
decoupled, because the resultant action in our theory 
does not contain the hybridization between bosons and 
fcrmions. The quasiparticlc Green's function is defined, 
as in the previous sections, as 



G { £(r) = -{Tf ia (r)fl(0)). (70) 

The bosonic part in Eq. (|69[) is given by 

4(r)% CT (0)) ~ 9LgL(Tlbl(T)-p i<T (r)} 

x [p] CT (r)-6 J -(r)]) ) (71) 

where we use vector notation as b i = (e|,dj), p\ a = 
{Pi/jiViS )■ Because we adopt the boson dynamics in 
which charge and spin bosons are decoupled, this bosonic 
part of the Green's function is rewritten as 

Tzl(r)z ja (0) 

'-91*92* ( T ibo -Po^iPoa -t>o 



+gtgL{T[b i (r)-pL}[Po*-b ] (r)} 
+glgl (T[b -p la {T)][pUr)-bl] 



where b = (e ,d ) and p 0a = (p 0a ,p ow ). If we retain 
only the first and second lines of the right hand side 
of Eq. (|72|) . the electron Green's function is reduced to 
that already obtained in Refl57l The contribution of the 
fourth line in the right hand side of Eq. (|72|) is small com- 
pared with these from other lines in Eq. ([72]) because it is 
the fourth order in the terms of the fluctuation denoted 
by tildes, and we ignore the fourth term. 



D. Superconductivities mediated by cofermions 
and charge bosons 

In this subsection, we examine how the present insight 
from the novel cofermion excitations offers the mecha- 
nism for the emergence of superconductivity. Here we ex- 
amine the singlet superconductivity induced by exchang- 
ing one charge boson (see FigE]) - 

■ (b) . I 



FIG. 6: Interactions between quasiparticles and cofermions 
mediated by exchanging one charge boson. 



We start from the mean-field action for the singlet su- 
perconductivity. Here the effective interaction of quasi- 
particles and cofermions is induced by exchanging one 
charge boson. The mean-field action Ssc is defined as 



4< 



S sc = T, *t(fe») [G{k,ie n )]- 1 * fc (fe n ), (73) 



9lo9lo (T[bl(r) • p i(T (r)][pt ( T ) . bJr)}) , (72) where #£(ie„) and $ fc (fe„) are defined as 



ft 



$ fc (i£„) 



f kt (is„), ^(ie n ) + ic^iien), f_ H (-ie n ),i>i ki (-ie n ) + xL H (-i£ n ) 



The inverse of the matrix form of the Green's function is given by 

is n -t k A k M k A fc 

A fc ikiEn-OL k Afc 

M k A fc iE n + Z k -Afc 
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(74) 
(75) 
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where = Co e k — A 4 , and M k and A k are superconducting order parameters as is explicitly given in Eqs. (|89|) and (|90[) 
Here the superconducting order parameters are given as 

T 



M k (ie n ) = — ^ t k t k+Q (b* Q (iui 7n )b Q {iui m ) 



and 



iuj m ,Q 

ep + da 
V2 



T 

A k (ie„) = "aT X] tktk+Q (b f Q (iuj m )b Q (iuj m ) 



[Q(k + Q,ie n + iu> m )] 13 + e ° ~t^° [<7(fc + Q,ie n + iw m )] 14 



(77) 



AT, 



[<?(fc + Q, ie„ + iw m )] 14 + e °^° + Q,te n + ^m)] 13 



(78) 



where [£(fc + Q, ie n ],y is the (i, j')-th element of the matrix Q. We should note that there are quasiparticle-quasiparticle 
and cofermion-quasiparticlc singlet pairings, which are described by following anomalous Green's functions: 



[Q(k,ie n )] 13 = 
[9(k,ie n )] 14 = 



2« fc A fe A fc - M k 


(tenf - a k 


(is n f - X 2 k+ 




fen) 2 - K- 



1 -M k A k (ie n -a k )+A k 


A fe + A l - ( i£ n - a k ) {ie n + &) 




(tenf ~ A2 + 




(tenf - \\_ 





where 



A fe ± 



X k ± y/X'i - AY k 

71 : 



Xu = 2jE+# + 2A£ + 2A£ + 2M£ 

2 



Y k = (a k M k + 2AA fc ) + [a*& + (k\ - A 
For simplicity, we introduce the cut off frequency 



uj c =max Cot, -v/ A fc + M fc 



(79) 
(80) 

(81) 

(82) 
(83) 

(84) 



to the bosonic propagators for charge fluctuations b\ = ej + d i and apply a quasi-static approximation to solve the 
gap equations as 

2 

M k {ie n ) x 



e + d 



V2 



[Q(k + Q,«£„ + iw m )]i3 + 6 ° "t-^ + Q,te n + iu m )] u 



+ 



— ^t k t k + Q (b f Q (uj c )b Q (u} c )^Tj2 
e -Mo 



V2 

—i \ 2 

ep + dp 

~7T 



(85) 



[Q(k + Q, ie„ + iu) m )] 



13 
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and 



A k (ie n ) x 



1 1 



V2 



(86) 



(87) 



[<?(fc + Q,i£ n + «*0] 14 + e °^° {&( k + Q> i£ n + iu> m )] 13 



15 



Then we obtain the superconducting states induced with electron-electron and cofermion-electron pairs formed by 
exchanging one charge boson. 



We solve the gap equations Eqs. ([55|) and ([55} with 
an assumption for the symmetry of the superconduct- 
ing order parameters. Here we assume a simple d x i_ y i- 
symmetry, which has experimentally been suggested for 
the hole-doped cuprates, for the singlet electron-electron 
and cofermion-electron pairs as 

Md 2 _ 2 (ie n ) 
M k (ie n ) = ^—^ (cosk x - cos fey) , (89) 

~ Ad 2 _ 2 (l£ n ) 

A fe (ie„) = — - — v - (cosfe^ - cosfey) . (90) 

As a consequence of the quasi-static approximation, the 
gap amplitudes Md a2 _ 2 (ie n ) and Ad x2 _ y2 (ie n ) become 
constant independently of the Matsubara frequency. Of 
course, we expect that the amplitudes of Md 2 2 (*£n) 

and Ad i2 _ y2 (is n ) will damp for lo c <C \ie n \. 

Here we note how the cofermion-quasiparticle pairs, 
namely, [£7] 14 , enhance the quasiparticle-quasiparticlc 
pairing. First, it is clearly seen in Eq. ([77j) or (|55|) . 
that the cofermion-quasiparticle pairs contribute to the 
superconducting order parameters of the quasiparticlc- 
quasiparticle channel, Mk(ie n ). This contribution 
becomes constructive one, only when the phases 
of the cofermion-quasiparticle and the quasiparticlc- 
quasiparticle pairs are the same, namely, Mk ■ Afc > 0. 
In this case, the interaction illustrated in FigJS^b) be- 
haves as an attractive interaction between the cofermion- 
quasiparticle and the quasiparticle-quasiparticle pairs. 
Then the formation of the cofermion-quasiparticle pairs 
enhances the quasiparticle-quasiparticle pairing through 
this attraction, and, indeed, enhances the pairing in our 
self-consistent solution of Eqs. (|8"o]) and (|88[) . 



IV. RESULTS 

In this section, we show how our theory predicts phys- 
ical properties of our interest. We estimate the quasi- 
particle and electron Green's function obtained in the 
above section, Scc lIII Cl We show the results for the 
Hubbard model defined in Eq.([T]) on a square lattice at 
U = 12t and t' = 0.25i to get insight into the cupratc 
superconductors. Wc restrict the mean-field solutions to 
the homogeneous and paramagnetic ones. All the calcu- 
lations are done at zero temperature. At this parameter 
and within the present calculation, the amplitude of the 
Hubbard gap SU at half filling n = 1 is estimated to be 
3.6t. 

First, we give the spectral functions calculated from 
the electron Green's function given in Eq. (|69p . and show 
the global structure of the spectra obtained from the elec- 
tron Green's function. Wc define the retarded Green's 



function at the frequency u> as 

G a (k,u)=g„(k,ie n ->u + S), (91) 
where 5 — > +0. Then, the spectral function is given by 

A(k,u) = lm[G a (k,u)]. (92) 




(0,0) (n,n) k (71,0) (0,0) 



FIG. 7: (color online) Spectral function A(k,ui) along lines 
running from (0,0) to (n,n), from to {it, 0), and from 

(tt,0) to (0,0). Here we use a finite broadening factor S — 
0.05t. 

In FigJTJ we show the result for x = 0.05. There are 
two main features. The first is the coherent band seen 
around the Fermi level, i.e., u = 0. The second is the in- 
coherent spectrum induced by the charge and spin boson 
dynamics. As is already discussed in RefloTl the remnant 
of the upper Hubbard band is seen above u ~ 6t, while 
the lower Hubbard band is seen between the coherent 
band and u> ~ —8t. The gap seen in the spectral func- 
tion at 2t < u> < 6t is nothing but the remnant of the 
Hubbard gap. Here an additional incoherent band is seen 
just above the Fermi level up to w ~ 2t. This additional 
incoherent band originates from the dynamics of the spin 
bosons. The spectral function shown in Fig(7] is qualita- 
tively consistent with the results of the quantum Monte 
Carlo simulations done by Preuss et al££, although their 
simulations were done for t'/t = 0, U/t = 8 and T ^ 0. 
We note that the dispersion of the coherent band is given 
by poles of the quasiparticle Green's function. Below, we 
concentrate on the low-energy excitations within the co- 
herent band. Therefore, we focus on the quasiparticle 
Green's functions instead of electrons Green's functions. 

In FigEl we show the LUSW defined by (1 + x)Co as 
a function of x. The quick increase of the LUSW larger 
than 2x is clearly seen in FigJ5J We refer to the data 
on the effective electron number N c g estimated from the 
optical conductivity measuremen t 37 ! 67 , which shows nice 
agreement with the present LUSW at the small doping 
as is expected. For larger doping the LUSW and N e g be- 
come trivially different because A c ff is also proportional 
to the population of the occupied states and decreases 
when the number of occupied states becomes small, while 
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FIG. 8: (color online) Bold (red) solid curve shows 
doping dependence of low-energy spectral weight. Thin 
(black) solid line represents 2x and thin (black) dashed 
line stand for x. O pen circles and closed squares illus- 
trate Nes given in Refl37l for La2- a; Sr :E Cu04, and Refl67l for 
Ca2- a: Na a ;Cu02Cl2 Respectively. 

the LUSW expresses only the population of unoccupied 
states and should always be larger than 2x. LUSW and 
N e g become identical in the small doping asymptotically. 



is 



where d k = § (( e k - fi± a k /j k ). 

Green's function given in Eq.(|93|) shows the divergence 
of the quasiparticle self-energy given by A|/ (7fcW — a k ), 
at 7^0; — a k — 0. In other words, the zero surface defined 
by 7fcW = a k emerges. Then, the zero surface splits the 
band dispersion defined by uj = Qo a e k — \i into two bands 

asw = d^±^J {d^) 2 + A^/7fc, as is depicted in FiglHl For 
small doping such as x = 0.05, our theory predicts that 
the reconstructed Fermi surface becomes a small pocket, 
as is seen in the right panel of Fig GO 

To show changes in the Fermi-surface topology with 
increasing doping, the single particle spectral function 
at uj = is given for the hole-doping rate, x = 
0.05, 0.15, 0.20 in FigHD{a)-(c), with 8 = 0.05*. The 



The quick increase of LUSW indicates that the Mott gap 
collapses and interrupted by the quick emergence of the 
low-energy unoccupied states. This LUSW comes from 
the quasiparticle band hybridized with the cofermions. 
The quick increase of the LUSW is also related with a 
positive feedback of the doping, where the hole doping in- 
troduces additional screening of carriers leading to a fur- 
ther increase of unbound holon and doublon. Although 
the Mott transition is continuous in the low-energy limit 
at the Fermi level as the ground state properties, this 
positive feedback gives a character close to the first-order 
transition in the energy scale of the LUSW, namely in the 
low-energy excitation spectra. It may be related to the 
tendency for the phase separation universally suggested 
in the doped Mott insulators ( including the cuprate su- 
perconductors ) which has been discussed from several 
different viewpoints^. 



A. Fermi-surface topology 



I 

topological transitions occur at x ~ 0.13 and x ~ 0.18, 
as is depicted in Fig lTUt e). In the region 0.13 < x, only 
small Fermi pockets exist, where the zero surface smears 
out the outer part of the pocket yielding an arc struc- 
ture as we see in FigfTOTa) (note that the damping is 
large near the zero surface because of an enhanced self- 
energy), consistently with the experimental signature in 
Fig. HOT d) as we discuss below. For 0.13 > x, a com- 
pletely different topology with large Fermi surfaces ap- 
pear, instead of Fermi pockets. For 0.13 < x < 0.18, a 
hole-like surface centered at (tt, 7t) (depicted by the solid 
curve in Fig fTOT b)) and, an electron- like one centered at 
(7T, 7r) (depicted by the dashed curve in FigflUTb)) coex- 
ist. These two surfaces together figure out an enclosed 



Now we show how the reconstruction of the Fermi surface occurs in our theory. The quasiparticle Green's function 

Gp>(k,u) = g<P(k,ie n ->u + iS), 



given as 



Gi f) (k,0j) = OJ + id-Coatk + V 




A 2 fe 



7 fc (cj + iS) - a k 



(93) 



^(d-y + Ai/ lk 



u + i5-d+-^{d- k Y + tf k h k 

1 

ui+is-di+^d^r+Ai/^ 



(94) 
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FIG. 9: (color online) Left panel shows band dispersions and 
zero surface calculated for x — 0.05 along lines running from 
(7r,7r) to (0,0), from (0,0) to (tt,0), and from (tt, 0) to (0,0). 
Right panel shows bare and reconstructed Fermi surface, and 
zero surface at cj = 0. Thin solid (black) curve gives the bare 
band dispersion to = C,o<r^k — Mi an d thin dashed (black) curve 
represents the zero surface 7^,0; = a^. Thick (blue and red) 
solid curves stand for reconstructed bands. 
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faces with the increasing electron concentrations. When 
the Fermi level is shifted by doping and touches the 
saddle points at (±7r, 0) and (0, ±7r), such topological 
changes occur. 

However, our solution shows another nontrivial topo- 
logical transition at x ~ 0.13. The topology in the phase 
x < 0.13 is highly nontrivial, and is not adiabatically 
connected with the conventional Fermi liquid. Such a 
topological change never occurs without the zeros of the 
Green's function. 

For comparison with experiments, we refer to the 
ARPES spectrum of Laa^Sr^CuO^ in FigOJId). Arc- 
like Fermi surfaces observed for x = 0.03 and x = 0.07 
show overall consistency with our result for x = 0.05. 
On the other hand, so-called nodal metallic behaviors, 
i.e., large amplitude of the spectral weight allowed only 
around the nodal direction (line running from (0, 0) to 
(tt, 7r)) are not seen in our results. Instead, in our the- 
ory, the weight is rather larger near the endpoint of the 
arc than the nodal point as we see in FigfTOTa). This 
subtle discrepancy is likely to originate from additional 
self-energy effects, which are not included in the present 
theory For example, fluctuations of d- wave superconduc- 
tivity and/or short range antifcrromagnctic fluctuations 
are possible candidates. However, we note that these 
fluctuations with finite correlation length by themselves 
never induce changes in the Fermi-surface topology. 



FIG. 10: (color online) (a),(b),(c) Single particle spectral 
function at w = 0. Solid and dashed lines illustrate the poles 
of quasiparticles. (d) Fermi surface observed by ARPES^. 
(e) Phase diagram of Fermi-surface topology in our theory. 
Topological phase boundaries exist at the doping x — 0.13 
and 0.18. 
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hole (electron unoccupied) strip between the solid and 
dashed curves. However, the electron-like surface (the 
dashesd curve) is hardly seen because the zero surface 
exists near this electron-like surface. On the other hand, 
for 0.18 > x, there exist an electron-like surface cen- 
tered at (0, 0) (depicted by the solid curve in FigfTUTc)). 
and an electron-like one centered at (tt, tt) (depicted by 
the dashed curve in Fig fTOT c)). The electron- like sur- 
face centered at (tt, tt) becomes less-visible for the used 
broadening factor S = 0.05t than the electron-like surface 
centered at (7r,7r) for x < 0.18. 

The quantum transition at x ~ 0.18 is a trivial one 
expected from the single particle picture. The topology 
of the Fermi surface at this quantum transition point 
changes from a hole-like surface centered at (tt, tt) de- 
picted for x = 0.15 in FigJTU] by the solid curve, to 
an electron-like surface centered at (0, 0) depicted for 
x = 0.20 in FigJTU] by the solid curve. This topology 
change is essentially understood from the non-interacting 
picture, where the band dispersion = — 2t(cosk x + 
cos k y ) + At' cos k x cos k y with t'/i = 0.25 shows the tran- 
sition from the electron-like to the hole-like Fermi sur- 



FIG. 11: (color online): Single particle spectral function near 
antinodal point for x = 0.10. Thin solid and dashed lines 
illustrate poles of quasiparticles. White bold dashed line 
stands for ARPES spectrum observed in L^-^Sr^CuO^ for 
x = 0.1022. 

When the large Fermi surface appears for x > 0.13, 
the Fermi surface always extends across the Brillouin 
zone boundary or the lines which connect the T point, 
(0,0), and the antinodal points (±7r, 0) or (0,±7t). How- 
ever, the single-particle excitation from the Fermi level 
now has a gap due to the hybridization gap Afc in ([93)1 in 
the quantum phase characterized by the nontrivial Fermi 
surface topology for x < 0.13. In FigfTTl the single par- 
ticle spectral function is depicted along the symmetry 
lines in the Brillouin zone. A gap measured from the 
Fermi level emerges, which corresponds to the pseudogap 
in the ARPES measurements. Hereafter, we define the 
amplitude of the pseudogap in the single particle spec- 
trum, A, as the gap amplitude between the Fermi level 
fi and the maximum of the single particle dispersion be- 
low (i along the line running from (0, 0) to (tt, 0) and 
the line running from (tt, 0) to (tt, tt). We also repro- 
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duce the ARPES dispersion observed in I^-xSr^CuCU 
for x = 0.10^ for comparison in FigfTT] (indicated by 
bold dashed curve). Here we employ a widely accepted 
parameter t =400 mcV. Indeed the dispersion has an ex- 
cellent similarity The quasiparticle dispersion does not 
touch the Fermi level along the Brillouin zone boundary 
here the line connecting (jt, 0) and (rr, tt), in both of our 
result and the ARPES measurement. In other words, the 
pseudogap amplitude A is finite for both cases. However, 
the quasiparticle state at (tt, 0) in the ARPES data has 
lower energy than our result and the discrepancy is up to 
20 meV. We expect some additional factors to push the 
dispersion down around the antinodal points (±7r, 0) and 
(0, ±7r). As is discussed above, fluctuations of <i-wave 
superconductivity and/or short range antiferromagnetic 
fluctuations are possible candidates of the origin of this 
discrepancy. 



to (tt, tt)) than around the antinodal points (±7r,0) and 

(0,±7T). 




(71,0) 



FIG. 12: Band dispersion of quasiparticles for x — 0.05. A 
full gap structure is clear. 



B. Pseudogap 

The pseudogap in single particle excitations defined in 
the above section characterizes a quantum phase with 
the nontrivial Fermi surface, because a nonzero pseudo- 
gap induces the breakdown of the trivial large Fermi sur- 
face. In our theory the pseudogap is induced by the 
hybridization gap between the cofermions and quasipar- 
ticles. Therefore, the amplitude of the pseudogap is 
roughly determined by the hybridization gap A& / ypfk. 
The explicit evaluation of the amplitude of is done 
by using Eq. (|B3|) in AppendixIBl in the lowest order of 
tij. The amplitude of A^ is roughly estimated as t 2 /W 
multiplied by a numerical factor, where W is the band 
width of original electrons. Here we note that the denom- 
inators in Eq. (|B3[) include the dispersion of the charge 
and spin bosons. The band width of the charge and spin 
bosons arc roughly equal to the band width of original 
electrons, W ~ 8t, and the band width of the spin wave, 
J ~ At 2 /U, respectively. For U = 12t, the band width of 
the charge bosons dominates that of the spin bonosons, 
namely, W > J. Therefore, the denominators appearing 
in Eq. (|B3|) are roughly scaled by W . On the other hand, 
the amplitude of 7/. is roughly estimated as t 2 /W 2 mul- 
tiplied by a numerical factor, through Eq. (|A17p . There- 
fore, the amplitude of the hybridization gap Ak/yfjk is 

roughly estimated as (t 2 /W)/ ' \Jt 2 /W 2 = t multiplied by 
a numerical factor. The present result for the energy 
scale of the pseudogap is consistent with previous stud- 
ies based on the cluster perturbation theory^ and, later, 
cellular DMFTZS. 

We discussed the pseudogap A along the definition of- 
ten used in ARPES measurements, in IIV Al Here we 
discuss the behaviors of the hybridization gap between 
two bands given in Eq. (|M]) and illustrate in FigGO First, 
the hybridization gap Afc / jk m our theory does not have 
nodes, and is s-wave-like in the momentum space, al- 
though, as is illustrated in FigfT2l the gap seems to be 
smaller in the nodal direction (the line running from (0, 0) 
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FIG. 13: Density of states (DOS) of electrons given as a func- 
tion of hole doping rate x. The thin dashed curve stands for 
the DOS of the noninteracting case. Closed (blue) circles show 
7 extrapolated from experimental data of La2-^Sr a: Cu04 in 
RefFrH. Closed (red) squares stand for experimental data of 
La2- a; Sr a ;Cu04 in Refl72l. Inset shows an example for extrap- 
olation used here, see the text below Eq. (|95p for details. 

Below the transition point x ~ 0.13, the density of 
states (DOS) of the original electrons cj^ at the Fermi 
level, /Of, defined as 



Pf 



-lmG ftJ {k,uj = 0) 



(95) 



is clearly suppressed, as is illustrated in Fig|T3l 
We compare p-p with the specific heat measured for 
La2- a; Sr a ;Cu0 4 71 i 72 . The specific heat coefficient 7 is 
expected to be proportional to the DOS. In the Som- 
merfeld's free electron model, the specific heat coefficient 
7 at T = is given by 7 = TT 2 pp/3. Here we use this re- 
lation and t =400 meV to calculate the expected values 
for 7 with our theory. 
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We note that 7 in Fig[13] indicated by closed (blue) 
circles is obtained through a linear extrapolation from 
low-temperature data in Refl7ll. and shows the possible 
lower limit for 7. From the data in Refl7ll. 7 in the un- 
derdoped region shows three characteristic temperature 
ranges. For example, the data for x = 0.08 (see the inset 
of FigH31) show (1) for 100 K< T < 200 K, a linear- 
temperature dependence, (2) for 50 K < T < 100 K, 
another linear-temperature dependence but with differ- 
ent slope from (1), and (3) for T < 50 K, the BCS-like 
jump. We linearly extrapolate 7 for 50 K < T < 100 K 
and obtain 7 at T = 0. On the other hand, data for 7 
indicated by closed (red) squares^ in Fig[13] may set an 
upper limit for 7, because, to exclude the influences from 
superconductivities, the authors of Refl72l used extrapo- 
lation from Zn-doped samples, which tend to show larger 
7 than samples without Zn-doping. 

Our result for 7 is consistent with the experiments, 
although we do not choose parameters specific to 
La2- 2; Sr 2 ;Cu04. On the other hand, the extrapolated 7 
(closed (blue) circles in Fig fTo) indicates further reduc- 
tion of the density of states at the Fermi level. It sug- 
gests possible roles of fluctuations from the <i-wave super- 
conductivity and/or antiferromagnetism in lower energy 
scale, which is left for future studies beyond the scope of 
this paper. 



C. Superconductivities 
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FIG. 14: Doping dependence of excitation gaps. The thin 
(black) solid line stands the pseudogap amplitude in the nor- 
mal state. This amplitude shows the contribution that has 
nothing to do with the superconductivity but rather to do 
with the precursor of the Mott gap. The bold (black) dashed 
line stands the minimum of the single particle excitation gap 
in the superconducting state along the Brillouin zone bound- 
ary or the diagonal line connecting the P point and the antin- 
odal points. The bold (black) solid line shows a refined es- 
timate of the single particle excitation gap in the supercon- 
ducting state obtained from the pair potential mediated by 
exchanging one charge boson, where the potential is multi- 
plied by a scale factor (1 — £0)- The ^-dependence of ex- 
citation gaps observed by ARPES for Bi2Sr2CaCu2 0s+<j^ 
are indicated by the thick (blue) dotted line. The open 
square stands for the superconducting gap of optimally doped 
La2-iSr x Cu04 estimated by ARPES measurements^. The 
thin (black) dashed curve stands for the density of supercon- 
ducting electrons n s . 



In studies on high-T c superconducting cuprates, the 
mechanism of the superconductivity is of course the cen- 
tral issued. However, the consensus on the mechanism 
has not been reached after more than twenty years of the 
discovery. Here we show how the novel cofermion mecha- 
nism for superconductivities developed in Scc. lIIID works 
and reproduces experimental observations on the high-T c 
superconducting cuprates. 

We examine the doping dependence of the single par- 
ticle excitation gap in the superconducting state based 
on the formalism developed in Sec. MID . The result is 
shown in Fig03] (bold dashed line). We adopt, as a defi- 
nition of the amplitude of the gap in the superconducting 
state, the minimum of the single particle excitation gap 
along the Brillouin zone boundary and the line connect- 
ing the r point (0,0) and the antinodal points (±7r, 0) 
and (0, ±7r), which is the same definition for the pseu- 
dogap in the normal state. This definition of the am- 
plitude of the gap is consistent with the data-analysis of 
the ARPES measurements. As is shown in FigQJ] (bold 
dashed line), the single particle excitation gap is still fi- 
nite around x = 0.3 for the <i-wave superconducting state 
in the present treatment, which is quantitatively incon- 
sistent with the experimental fact that the superconduc- 
tivity disappears for x > 0.26 in the high-T c cuprates. 
To highlight a possible origin of this inconsistency, we 
also show a result for the single particle excitation gap in 
FigHH (thick solid line) calculated by rescaling the quasi- 



static boson propagator as 

fqt^^k)) -> (1 - Co) (^j(w c )bg(wc)) • 

By the rescaling, the single particle excitation gap is 
clearly suppressed around x = 0.3. The reason why we 
consider such a rescaling is the following: In our treat- 
ment, the charge fluctuations, especially holon fluctua- 
tions remain finite even in the dilute limit x — ^ 1, where 
holons fill all the sites and should be localized as hard 
core bosons. Holons are completely localized and form 
the "Mott insulating" state for x = 1. The fluidity of 
the holons for x < 1 originates from finite electron den- 
sity, 1 — x ^ 0, because the finite electron density intro- 
duces "vacancies" in the holon's Mott insulating state, 
and causes the fluidity of the holons. On the other hand, 
doublons vanish at the limit x — > 1. Near this limit, the 
density of the doublons are scaled by (1 — x) 2 . There- 
fore, for large doping x, the weight of fluctuating charge 
bosons should vanish as a function of 1 — x. We note that, 
for x ~ 1, 1 — Co — 1 — x. On the other hand, for x ~ 0, 
the charge fluctuating bosons have the weight 1 — Co, be- 
cause the incoherent band induced by fluctuating bosons 
share the total spectral weight with the coherent band, 
which has the weight Co- Therefore for all the doping 
range, the rescaling of the boson propagators is naturally 
interpolated by 1 — Co- 
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Then wc compare our result with experimentally ob- 
served gap A. The result of A without rescaling of the 
propagators given by the thick dashed curve in FigfT4l 
shows a rough qualitative agreement with the ARPES 
result obtained from the single particle excitation gap 
at low temperatures (T=14 K) for Bi2Sr2CaCu20s+^ 
(dotted (blue) line in Fig fTi)) . However, our result with- 
out the rescaling overestimates A around x = 0.3, where 
the superconducting phase is not observed in any of the 
cupratcs. On the other hand, the single particle excita- 
tion gap calculated by using the rescaled boson propa- 
gator shows a suppression of A around x = 0.3. The 
single particle excitation gap A for x ~ 0.1 is consis- 
tent with the superconducting gap of optimally doped 
La2- 2; Sr 2 ;Cu04 estimated by ARPES measurements^ 
(open square in Fig|TJ|). On the other hand, the rescaled 
boson propagator appears to be still overestimated for 
the overdoped region, x ~ 0.3. The origin of such overes- 
timated fluctuating charge bosons in our theory may be 
ascribed to the approximation on the bosonic propaga- 
tors given in the Dyson equations shown in FigjSJ Life- 
time of the coherent propagation of the fluctuating charge 
bosons is not taken into account in our theory. There ex- 
ist two possible origins of the finite lifetime. The first 
is damping caused by single particle excitations, namely, 
the Landau damping. When the doping increases, the 
motion of the coherent carriers may disturb the propa- 
gation of the fluctuating charge bosons seriously. The 
second is the separation of charge and spin fluctuations, 
which becomes a good approximation for x ~ 0. Such 
separation of charge and spin fluctuations may become 
worse for large x, because the motion of the coherent 
carriers mixes the charge and spin modes. The mixing 
of charge and spin fluctuations also disturb the propa- 
gation of the fluctuating charge bosons and give a finite 
lifetime. Therefore, these factors neglected in our theory 
will suppress the charge fluctuations and consequently, 
the superconductivities. 

The single particle gap A decreases monotonically, as 
the doping x increases. It is consistent with ARPES 
measurements. On the other hand, the critical temper- 
ature of the superconductivity, T c , shows a peak around 
the x ~ 0.15 in La2_ a; Sr 2 ;Cu04, and exhibits a dome- 
like structure as a function of x. To estimate the ten- 
dency of T c from our results for T = 0, we show the x- 
dependences of the density of superconducting electrons 
n a in FigQjJthin dashed curve) defined as 



Cojt- E ( cos k x - cos k y ) [g(k, ie n )' 
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(96) 



The pairing mechanism proposed here includes the 
contribution of the cofermion-quasiparticle pairs in ad- 
dition to the quasiparticle-quasiparticle pairs. This con- 
tribution of the cofermion-quasiparticle pairs itself is a 
novel perspective introduced in this paper. Further- 
more, it offers an insight into one of the most inter- 
esting issues in physics of the cuprate superconductors, 



namely, relationship between the pseudogap formation 
and the high-T c superconductivity. The cofermions in- 
duce the hybridization gap around the Fermi level, and 
as a result, the pseudogap. The pseudogap formation 
itself reduces the DOS around the Fermi level and is 
harmful for the high-T c superconductivity. However, 
the cofermions support the high-T c superconductivity 
through the cofermion-quasiparticle pairing, simultane- 
ously. Therefore the pseudogap structure does not nec- 
essarily destroy the superconducting pairing. The pseu- 
dogap formation is the other side of the coin of the 
cofcrmion pairing contributing to the superconductiv- 
ity. This dual character offers an in sight into the re- 
cent controversy called dichotomy observed in ARPES 
measurement a 38 ' 39 : There exist two major pictures for 
the relationship between the pseudogap and high-T c su- 
perconductivity. One tells us that the pseudogap is a 
precursor of the Mott gap£^, and coexists with the su- 
perconducting gap. The other claims that preformed 
pairs of strong coupling superconductivities induce the 
pseudogap^. The present theory offers a new alternate 
route of understanding that reconciles the dichotomy. 

By exchanging one charge boson, attractive interac- 
tions are induced for forward scattering channels. Even 
though there is the strong on-site repulsion U, such at- 
tractive interactions may favor anisotropic superconduc- 
tivities not only the d a: 2_ y 2-wave superconductivity but 
also, for example, d xy -wave and extended s-wave super- 
conductivities. In the present theory, however, a sim- 
ple extended s-wave superconducting parameter propor- 
tional to cos k x + cos k y does not develop. 

D. Metal-insulator transitions 

We propose a scenario for the filling control metal- 
Mott insulator transitions based on our theory. Before 
going into our scenario, we note that there is a difficulty 
in our theory in the small doping limit x — > +0. In 
this limit, we have Fermi pockets with small but finite 
volume depending on parameters such as t'/t and U/t. 
It is in contrast to the zero-doping limit of the ordered 
insulators. In our theory, the weights of the two bands 
split by the zero surface illustrated in Fig|5] are not the 
same. From Eq. ([M)l . these weights are given as 
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In the case of the antiferromagnetic ordered state, the 
counterparts of these weights are given by 
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from Ea.(|49]). and are the same, because of 
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The weight compensation between these two bands split 
by the zero surface is needed to reproduce ordered insu- 
lating phases. On the other hand, our theory does not 
offer this weight compensation. Therefore, our theory 
does not necessarily predict that the small hole pockets 
seen in Fig lTOt a) shrink to points at the small doping 
limit x — > +0. 

In addition, in the limit x — > +0, the imaginary part 
of the quasiparticlc and cofermion self-energy, which is 
not taken into account in the present theory, may affect 
the weight compensation. Higher order corrections to the 
cofermion dynamics will also affect the structure of the 
zero surface in this limit. 

Despite these limitations, our result still suggest a 
possible scenario for the filling control metal-Mott in- 
sulator transitions. If the small Fermi pockets seen in 
Fig fTor a) shrink to the Fermi points, the filling control 
metal-insulator transition occurs as a topological transi- 
tion. On the other hand, if the small Fermi pockets do 
not shrink to the Fermi points, the quasiparticle weight 
of the Fermi pockets becomes zero, in spite of the finite 
volume of pockets. Then the insulating state appears, 
because of the fading out of the coherent band. In both 
of the two cases, a small but finite DOS is expected in 
the limit x — > +0. It should be noted that, anyway, a 
nontrivial topological transition is expected to occur at 
a finite doping, before the system becomes insulating. 

Here we note that the order of the three quantum phase 
transitions, namely, the Mott transition at x = 0, and 
the two topological transitions at x ~ 0.13 and x ~ 0.18 
are all continuous within the present calculated results, 
although in general, the first order transition is not a pri- 
ori excluded. If a topological quantum phase transition 
takes place as a continuous one at T = 0, it becomes just 
a crossover at nonzero temperature, because the Fermi 
surface is in the strict sense not well defined at T > and 
the toplogical character of the Fermi surface loses well- 
defined meaning at T > any more. It has been stressed 
before that the sharp change between the overdoped and 
underdoped regions is associated with some kind of quan- 
tum criticality triggered by a quantum phase transition 
of the symmetry breaking order. In this conventional 
quantum criticality, the phase transition usually extends 
to nonzero temperature as the border of symmetry bro- 
ken phase. However, the present quantum criticality is 
completely different from this category, where not the 
symmetry breaking but the topology change drives the 
transition's. 



V. SUMMARY AND DISCUSSIONS 

In this paper, we have proposed a microscopic mecha- 
nism for the Fermi-surface reconstruction and non- Fermi- 
liquid behaviors emerging in proximity to the Mott in- 
sulators, motivated by puzzles experimentally observed 
in the underdoped cupratcs. We construct our theory 
starting from one of the simplest theory for correlated 



metals, namely, the slave-boson mean-field theory for 
the Hubbard model by Kotliar and Ruckenstein (KR). 
Then a special emphasis is placed on the role of extra 
charge dynamics on low-energy spectra. We find that 
hidden cofermionic particles constructed from composite 
fermions play a crucial role in addition to the quasiparti- 
cles described by the previous KR theory. The additional 
cofcrmions called holo- electrons and doublo-holes repre- 
sent substantial low-energy part of the charge dynamics. 

Thus introduced new cofermions hybridize with the 
quasiparticles and cause a hybridization gap represented 
by zeros of the quasiparticle Green functions. As a re- 
sult, although the explicit symmetry breaking is absent, 
the gap identified as the pseudogap of the cuprates nat- 
urally comes out. Although the origin of the gap has an 
apparent similarity to the case of the symmetry-broken 
phases such as commensurate antiferromagnetic metals, 
in the sense that both of them can be ascribed to a hy- 
bridization gap, the mechanism of the hybridization and 
the partner of the hybridization are very different from 
the simple symmetry breaking. The origin of the gap is 
also different from the Mott gap itself. The large Mott 
gap in the insulating phase is quickly replaced with the 
present much smaller hybridization gap upon doping, be- 
cause the LUSW above the hybridization gap emerges. In 
the result of the cellular dynamical mean-field theory^, 
this hybridization gap coexists with the clear remnant of 
Mott gap structure in a much larger energy scale. This 
indicates that the two gaps are separated phenomena. 
Thus the pseudogap or the hybridization gap in the un- 
derdoped region is not really the precursor of the Mott 
gap. 

When the cofermion strongly hybridizes with the 
quasiparticles and the resultant hybridization gap splits 
the quasiparticle band around the Fermi level, Fermi- 
surface reconstruction occurs in a natural way. The hy- 
bridization gap indeed leads to a break-up of the Fermi 
surface into small hole pockets around the nodal direc- 
tion, and a phase topologically different from the predic- 
tion of the weak coupling picture emerges. This phase re- 
alized in the underdoped region is separated by a topolog- 
ical phase transition from the normal Fermi liquid phase 
in the overdoped region. Our result clearly shows that 
an unconventional non-Fermi liquid phase emerges in the 
underdoped region. Such topological changes and an 
emergence of a new phase are consistent with experimen- 
tally observed "Fermi-arc" formation in hole-undcrdopcd 
cuprates. 

The topological quantum phase transition between the 
conventional Fermi liquid and the non Fermi liquid is a 
continuous one at T — within the present approxima- 
tion. This implies that it transforms to a crossover at 
nonzero temperature, because the continuous topologi- 
cal transition is well defined only at zero temperature. 
In principle, it does not exclude a possibility of the first 
order transition of the topological change, which extends 
to nonzero temperatures with the critical end point at a 
nonzero temperature^. The present results do not sup- 
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port the existence of such a first order transition in agree- 
ment with the absence of the clear indication of the first 
order transition in the experiments. 

Our formulation offers a new concept for the "pseudo- 
gap" phenomena observed in hole-underdoped cuprates. 
It proposes that the "pseudogap" in the single particle 
spectrum results from the hybridization gap between the 
cofermions and the quasiparticles. The part below the 
Fermi level has the maximum in the antinodal region 
and has a structure similar to the d x 2_ y 2 symmetry in 
agreement with the experimental observation. The main 
part of the gap lies, however, above the Fermi level. The 
total hybridization gap has the structure of the s-wave 
gap. 

Consequently, the characteristic energy scale of the 
pseudogap is identified as that of the hybridization gap 
amplitude and thus basically scaled by the bandwidth 
(kincic energy) of the bare electrons, t multiplied by a nu- 
merical factor. Therefore, the energy scale of the "pseu- 
dogap" has something to do neither with the Hubbard 
gap or the on-site Coulomb repulsion U nor with the su- 
perexchange interaction J ~ t 2 /U. 

We have also numerically studied the relevance of our 
theory in the realistic condition. For the on-site Coulomb 
repulsion U large enough to create the Hubbard gap 
?7 > lOt and small amount of hole doping x < 0.1, re- 
garded as a relevant parameter for the high-T c cuprates, 
the structure of the momentum dependence and the am- 
plitude of the pseudogap below the Fermi level have a 
quantitative consistency between the ARPES measure- 
ments for the hole-doped cuprates and the present re- 
sults. Doping dependence of the pseudogap amplitude is 
also consistent with the experimental observation. 

Reduction of the electronic density of states induced by 
the pseudogap formation gives another consistency with 
the experiments as observed as specific heat coefficients. 
The pseudogap formation in our theory also induces 
asymmetry of the DOS around the Fermi level, which 
naturally explains the asymmetric STM spectr a 81 ' 82 . 
The overall consistency supports that the cofermions, 
holo-electrons and dublo-holes, are indeed a relevant ob- 
ject to be considered in physics of the cuprate supercon- 
ductors. 

The present theory has further consequences and pre- 
dictions for experiments. The hybridization gap basically 
has a s-wave- like symmetry and a major part of the gap 
lies above the Fermi level for the model of the hole doped 
cuprates. Although the present resolution limit of the in- 
verse photoemission does not allow determination of the 
detailed structure of the unoccupied electronic states, we 
propose that our prediction of this s-wavc-likc gap struc- 
ture accompanied by the main part of the LUSW lying 
above the hybridization gap can be tested experimentally 
if some high-resolution measurements of the unoccupied 
spectra are provided. The optical conductivity tr(w), es- 
pecially mid-infrared peak and a long tail of <j(oS) ob- 
served in the underdoped cuprates indeed supports the 
existence of such LUSW— iSI. 



So far, direct measurement of the cofermions appears 
to be difficult, because the dispersion of the cofermions 
corresponds to the zeros of the quasiparticles, while the 
zeros are in general hard to detect. In addition, the 
cofermion does not have a charge and does not allow an 
electromagnetic detection. However, it contributes to the 
entropy and thermal transport. Since the electric conduc- 
tion is contributed only from the quasiparticlc, while the 
thermal transport can arise from the cofermion as well, 
we expect a serious breakdown of the Wiedcman-Franz 
(WF) law. The WF law predicts that the ratio of the 
thermal conductivity k to the electric conductivity a is 
equal to LqT, namely, Lq = n/Ta, where Lq is a uni- 
versal constant given by L = {ir 2 /3) ■ (ks/e) , namely, 
the Lorenz number. On the other hand, in our theory, 
the ratio L = n/Ta is predicted to be larger than Lq 
in proximity to Mott insulators, because the cofermions 
carry additional energy, and contribute to k. Indeed the 
breakdown of the Wiedeman- Franz law reported recently 
in hole-underdoped cuprates^ supports the present pre- 
diction. 

We propose a scenario for the filling control metal- 
Mott insulator transitions in two dimensions based on 
the newly introduced cofermions. We predict a pseu- 
dogap formation and consequently a non-trivial change 
in the Fermi-surface topology in the underdoped region. 
Then the criticality of the small Fermi pockets determines 
the nature of the Mott transition, if the transition is con- 
tinuous. Two scenarios remain possible: On the verge of 
the Mott transition, it may either shrink to points before 
vanishing by keeping the quasiparticle weight nonzero, 
or the quasiparticle weight decreases to vanish by keep- 
ing a finite volume of the Fermi surface. In the former 
case, the filling control metal-Mott insulator transitions 
occur as a topological one. The density of states remain 
nonzero and finite on the verge of the transition because 
of the two deminsionality. In the latter case, the effec- 
tive mass diverges. These two possibilities are essentially 
identical to the two types discussed in the literature^. 
Although the present numerical accuracy is not sufficient 
for determining the ultimate criticality, the shrinkage of 
the pocket and arc overall supports the former scenario, 
while the reduction of the DOS suggests only a part of 
the pocket contributes to low-energy excitations as the 
arc. This agrees with the experimental results. More de- 
tailed and accurate determination of the criticality is left 
for future studies. 

We also propose a novel mechanism for high- 
temperature superconductivity driven by the cofcrmion- 
quasiparticlc pairing. It offers a new insight into the 
relationship between the pseudogap formation and the 
high-T c superconductivity. The cofermions newly intro- 
duced in the present paper induce the hybridization gap 
around the Fermi level, and as a result, the pseudogap. 
This pseudogap formation itself reduces the DOS around 
the Fermi level and destroys superconductivity. On the 
other hand, the cofermions enhance the high-T c super- 
conductivity through the cofermion-quasiparticle pair- 
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ing, simultaneously. At the present level of approxi- 
mation, in the overdoped region, a process exchanging 
one charge boson overestimates the charge fluctuations 
and predicts much larger single particle gap in supercon- 
ducting phases, than experimentally observed gaps. It 
is left for future studies to correct such an overestimate 
of charge fluctuations. In addition to one charge boson 
exchange processes, cofermion polarization will help the 
pairing between quasiparticles, which is not taken into 
account as a higher order contribution in the present pa- 
per. However, the superconducting gap obtained from 
the present approximation already reproduces the dome- 
like structure as a function of the doping concentration 
with a right order of magnitude if we compare with the 
cuprates. 

Although we focus on the hole-doped cuprates in this 
paper, electron-underdoped cuprates are also interesting 
from the viewpoint of the present cofermion theory. In 
the electron-underdoped systems, our theory will predict 
the emergence of two electron pockets centered at (tt, 0) 
and (0, 7r) even in the absence of antifcrromagnctic long- 
range orders. Such a Fermi-surface topology may reduce 
instabilities towards antifcrromagnctic orders, and will 
make c? ;c 2_ y 2-wave superconductivities stable: The elec- 
tron pocket formation prevents the nesting of the Fermi 
surface, and allows an emergence of full-gapped super- 
conductors but with the (4,2 _ y 2 -wave symmetry, because 
the node lines run the momenta where the original elec- 
tron pockets are absent. 

Here, we make some remarks in regard to the relation- 
ship with other theoretical approaches. In contrast to 
our approach, the high-energy charge degrees of freedom 
to do with the UHB have often been integrated out to 
extract the low-energy physics from the Hubbard model. 
The t-J model is a typical effective model derived from 
such a treatment, and has been studied as a canonical 
model describing the doped Mott insulator. The gauge 
theory based on the t-J model has been intensively stud- 
ied to explain the low-energy physics of cuprate supercon- 
ductors as doped spin liquids^. Wen and Lee 8 -^ proposed 
that a hole-doped spin liquid state may exhibit hole pock- 
ets. According to this theory, the optical or direct gap of 
this phase is d-wavc-like: There is no optical gap along 
the nodal direction. Although the gapless excitation is 
present in the nodal direction by the electron-hole exci- 
tation through the hole pocket even in our theory, the 
gauge-theory scenario by Wen and Lee is in contrast to 
our result predicting that a relatively small but nonzero 
amplitude of the optical gap even along the nodal direc- 
tion is superimposed, namely, a s- wave-like optical gap 
should be visible between the two bands separated by the 
hybridization gap. The s- wave-like pseudogap is also sup- 
ported in the previous numerical studies by Stanescu and 
Kotliar— , and Sakai et al^, in support of the present 
theory. An earlier exact diagonalization study without 
any bias has indicated that even the doped t-J model 
has a similar s-wave-like pseudogap structure in agree- 
ment with our theory^. As we mentioned already, it is 



desired to test against the two contradicting predictions 
of our theory and the gauge theory by inverse photoemis- 
sion spectroscopies or by more sophisticated and high- 
resolution experimental tools to measure the unoccupied 
states in the future. 

Another crucial difference between the present theory 
and the gauge theory can be tested by the breakdown 
of the Wiedeman- Franz law. In the gauge theory, the 
breakdown is due to the contribution from the spinons^!. 
However, if the arc structure is observed, the d-wavc like 
gap needs to be developed by the flux fluctuations, which 
leads to the confinement of a holon and a spinon gener- 
ating a quasiparticle in the nordal direction. It should 
recover the Wiedeman-Franz law in this region of the arc 
formation. Therefore, at low temperatures in the under- 
doped region, the Wiedeman-Franz law should eventually 
be followed. On the contrary, in the present case of the 
cofcrmions, the contribution continues even at low tem- 
peratures and the breakdown of the Wiedeman- Franz law 
is robust. Recent experimental results appear to support 
our predictio n 83 ! 88 . 

Integrating out the "high-energy" charge degrees of 
freedom represented by the existence of doublon thor- 
oughly leads to the ignorance and overlook of impor- 
tant aspects of the low-energy spectrum of the Hubbard 
model. There exists a variety of experimental facts which 
are not accounted for by the t-J physics: An example 
is the estimate of the LUSW by using N e g(uj) observed 
in the optical conductivity measurement^!. Choy et al. 
have made a careful treatment to extract the low-energy 
effective action, with the failure of the t-J model kept in 
minds*. They claim the weight transfer among the co- 
herent band, the LHB, and the UHB described by hidden 
2e-bosons. In spite of the illuminating proposal, their ef- 
fective action could not be solved exactly. Approximate 
evaluations of their theory show a soft gap behaviors, 
and predict a semiconducting behaviors in the "pseudo- 
gap phase," where resistivity p is expected to diverge as 
p oc T _1 with lowering temperatures. In the present the- 
ory, the Fermi pocket appears in the "pseudogap phase" 
in the underdoped region, say, for x < 0.13. Therefore, 
we predict metallic conduction in the "pseudogap phase," 
in sharp contrast to the hidden 2e-boson theory. 

As the authors have already discussed in Refl27l. the 
formation of the cofcrmions shares profound similarity 
with that of the excitions in semiconductors. In the so- 
called d-p model for the cuprates, which contains both 
d-electrons on the copper sites and p-electrons on the 
oxygen sites, it has been claimed that the excitionic ef- 
fects due to d-p interactions strongly affects its excitation 
spectrum^. Relationship between exciton-like features 
of the cofermions and excitionic effects studied in the d-p 
model are also desired to be clarified in the future. 

Our theory has been constructed to account for the 
charge dynamics more seriously than the literature and 
proposed the topological changes of the Fermi surface 
in proximity to the Mott insulators. The reconstructed 
Fermi surface also offers an unexplored avenue at smaller 
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energy scale, if it is combined with other possible fluc- 
tuations such as spin and superconducting fluctuations. 
Clarifying possible emergence of antiferromagnetic orders 
in the small doping region x < 0.02 by using a unified 
scheme is left for future studies. It is also left for future 
studies how the topological changes affect various possi- 
ble symmetry breakings such as time reversal symmetry 
breakings, stripe formation, and incommensurate charge 
orders including phase separations. Experimentally sug- 
gested nodal metallic behaviors should also be examined 
in more detail. Our theory will give a new insight into 
emergence of these competing orders, the high-T c super- 
conductivity, and the anomalous metallic state. 

I 

Then the self-energy matrix in Eq.(|63|) is calculated by 
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Appendix A: Self-energy for the cofermions 

Here we calculate the self-energy for the cofermions in 
Dyson equations (see Fig[5]). 
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When we use the mean- field propagators for the quasiparticlcs, (f^ a f\ a \ wc obtain a simple analytic expression for 
the self-energy matrix as 
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Wf(Q) W* b (Q) 



>m ~ \£q\ 

Zf(ff) 



■n + \t Q \' 

Z ab (ff) 



iuj m -\A Q \ iw m + |Ag|' 
where the coefficient matrices, W± b and Z^, are given as 

1 f 1 A SX- ai^e Q / 1 0\ h^tQ ( I 



W l ± \Q) W?(Q) 
WlHQ) Wf{Q) 



±- 



—I 



2i c 



Z£(Q) Z¥(Q) \ _ SX + SU/2 ( 1 V 1 / 1 
Zf{Q) Zf{Q)) 2a Q [0 l) ± 2\0 -1 



21q V 1 °. 
c iTT e Q / 1 \ di^e Q ( 1 



2o c 



f 



2oq 



1 



(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



r 



The parameters used in the above equations are given 
as SX = A« - A< 2 ), SU = U - 2X^ 2 \ |e| = 



wXk,ie n ^QP(k,ie n 



SX - ai — e Q 



(A7) 
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°Q = 


/ 1 \di i 

V v 


Aq = 


su 




Aq = 


SU 





SU 
~2~ 



(A9) 
(A10) 

The coefficients, <xi, &i, ci, and di are determined as, 



zeoflo i \ e i a j I 






2 Po + (plaPja) " 




(A13) 


2 P0 + (PiaPja) "+ 




(A14) 



_ where we oniy take into account nearest-neighbor pairs 

^ +3 + (e«ej) + (S?e0 , (All) for 



J 



To derive Eq. (|A2[) . we also use following relations, 

1 1 



T E 



ie„ + iw m - £ iw m - Ai 



fl(0fl(-Ai)-fl(-Qfl(Ai) 
«£n - C + A 



(A15) 



and 



t3 EEE' 



1 



1 



ie„ + ibJra — iwt + iw„ — £ iw m — Ai iuig — A2 «w„ — A3 

g(Qg(-Ai)g(A 2 )g(-A 3 ) - g(-Qg(A 1 )0(-A 2 )^(A 3 ) 
i£„ - £ + Ai - A 2 + A 3 



(A16) 



Derivatives of E cd (fc) with respect to ie n give jk after taking a limit ie n — > 0, as is mentioned in Scc lIII Ci Here we 
give the derivatives as 



<9£ cd 



^ E E ^92*ft k+p t k+R (^ a] p J c p^(^ b R J d ^ Ra 



E 

Q 



0,6=1,2 



(l e fc+Q "Ml + l-M) (\ e k+Q - Ml + l A c 

i'k+Q P+R - t*)Z°*(Q)Wf(P)Wl d (R) 



7V3 E E C9i^2<x) 4 (ifc+tf) 



P,Q,R»*=l,2 



\ e k+Q P+R - Ml + |Aq| + IM + M 



(m - ^ + g-p+i?)^f(Q)^- b (^)w / + W 

e k+Q-P+R - Ml + |Aq| + Kp| + Kfll 



(A17) 



Appendix B: Hybridization 

The amplitudes of hybridization between quasiparticles and cofermions, 



A W ,A (X) 



are given as (A^,, Ai) = , A k )■ Its Fourier transformation is given as 
T 2 



A fc = a^EE E tk + R(9l3l) 2 (Pp-nP b p-n 

s P,R q o,6,c=l,2 



tfl ^ g+ fl . (B2) 



(Bl) 
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where we use vector notations, b = (6q,&q) = (eo,do), f3 t — (Pj,f3f) = (e i5 dj), and 4> i 
if we use the mean-field propagators for the quasiparticle, we obtain 



= ~J^2 J2 J2 t k+R(9ia9iaY ''tq+RPq 

s P,R,q a ' b ' c=1 ' 2 
' Zf{P - R)Wl b (P)Wl c (R) Z a }{P - R)Wf{P)W dc {R) 



\^P-R\ + \tp\ + \t R \ 



) = (P«tiJ4)- Here > 



(B3) 



where we use the following relation: Equation (|B3j) will help us to estimate the amplitude of 

the pseudogap in our theory. 

t 2 YY — 

z — 4 ^— ' icu£ — iuj n — Ai icu£ — A2 ico n — A3 

= fl(Ai)fl(-A 2 )fl(A 3 ) + 0(-Ai)g(A 2 )fl(-A 3 ) 
Ai — A2 + A3 
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